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ABSTRACT 


For shape predictions of structures under large geometrically nonlinear deformations, Curved 
Displacement Transfer Functions were formulated based on a curved displacement, traced by a material 
point from the undeformed position to deformed position. The embedded beam (depth-wise cross section 
of a structure along a surface strain-sensing line) was discretized into multiple small domains, with 
domain junctures matching the strain-sensing stations. Thus, the surface strain distribution could be 
described with a piecewise linear or a piecewise nonlinear function. The discretization approach enabled 
piecewise integrations of the embedded-beam curvature equations to yield the Curved Displacement 
Transfer Functions, expressed in terms of embedded beam geometrical parameters and surface strains. By 
entering the surface strain data into the Displacement Transfer Functions, deflections along each 
embedded beam can be calculated at multiple points for mapping the overall structural deformed shapes. 
Finite-element linear and nonlinear analyses of a tapered cantilever tubular beam were performed to 
generate linear and nonlinear surface strains and the associated deflections to be used for validation. The 
shape prediction accuracies were then determined by comparing the theoretical deflections with the finite- 
element-generated deflections. The results show that the newly developed Curved Displacement Transfer 
Functions are very accurate for shape predictions of structures under large geometrically nonlinear 
deformations. 


NOMENCLATURE 

c depth factor (vertical distance from neutral axis to bottom surface of uniform embedded 
beam), in 

c(s) = c(x) , depth factor, vertical distance from neutral axis to bottom surface of in- 
extensional (s = x) nonuniform embedded beam, in 

c = c(x;), depth factor at x = x,, in 

CG; =h,—c,, depth factor at x = x, for upper surface, in 

oF =c(x,), depth factor at X= x,, in 

Cp value of c, at embedded beam root x = x, =0, in 

C,, value of c, at embedded beam tip x= x, =/, in 


d math symbol—differential 

E Young’s modulus, Ib/in” 

h embedded beam depth at x = x,, in 
h 


0 value of h, at embedded beam root x= x, =0, in 
h,, value of h, at embedded beam tip x = x, =/, in 
l length of embedded beam, in 
n index for the last span-wise strain-sensing station (or number of domains) 
P applied load, Ib 
R(s) radius of curvature, in 
Ky curved axial coordinate along elastic curve of deformed embedded beam, in 
x,y Cartesian coordinates (x in beam axial direction y in lateral direction), in 
S; deformed curved axial coordinate of strain-sensing station at s,(= x,), in 
t tubular beam wall thickness, in 
U,; magnitude of displacement along x-axis of deformed material point at 5, , in 
Xx; x coordinate of the i-th strain-sensing station 


y, straight deflection at axial location x = x,, in 


LV poy value of y, calculated from Nastran linear analysis, in 

y(x) curved deflection (curved distance traced by a material point from its initial un- 
deformed position to its final deformed position), in 

y; = y(x,;) = y(s,), curved deflection at axial location x = x,, in 

Care value of y, calculated from Nastran nonlinear analysis, in 

y, beam-tip (i= 17) vertical deflection, in 

y, beam-tip (i= 7) curved deflection, in 

Al = (x, —x,_,) =(s,—5,_,)=1/n, domain length (strain-sensing stations separation distance), 
in 

€(s) surface strain at curved axial location s, in/in 

E(x) surface strain at x-location, in 

E; lower surface strain at i-th strain-sensing station, in/in 

E; upper surface strain at i-th strain-sensing station, in/in 

E, axial strain in s-direction, in/in 

A(s) beam slope angle in reference to s-system, rad or deg 

0(x) beam slope in reference to x-system, rad or deg 

6, = 0(s,) [= 0(x;)], slope angle at i-th strain-sensing station, rad or deg 

0, = 0(s,) [= O(x, )], slope angle at 0-th (beam root) strain-sensing station, rad or deg 

6. = 0(s,,) [= O(x,,)] , slope angle at n-th (beam tip) strain-sensing station, rad or deg 


slope angle of a straight line connecting origin and deformed material point s, on 


elastic curve of deformed embedded beam, rad or deg 
Poisson’s ratio 


um S 


=xX-—xX,_,, shifted axial coordinate, in 


INTRODUCTION 


After the Helios prototype with a wing span 247 ft (fig. 1) broke-up in mid-air at an altitude of 
2,800 ft under very large wing dihedral deformation with a wing tip deflection reaching 40 ft (fig. 2), 
there was a need to develop a new technology for monitoring the deformations of highly flexible aircraft 
for feedback control and flight safety. A new technology has been developed that can convert distributed 
surface strain data into structure deformed shapes. 

Strain sensors can only measure surface strains and not the structure deformed shape. However, after 
the development of the Ko Displacement Theory (refs. 1-8), a new innovative structure shape prediction 
technology, called Method for Real-Time Structure Shape-Sensing (U.S. Patent Number 7,520,176) 
(ref. 3), was created. This new technology uses the Displacement Transfer Functions to transform 
rectilinearly distributed surface strains into out-of-plane deflections for mapping overall structure 
deformed shapes for visual displays. For applications, the surface strains can be obtained from 
conventional strain gauges, wireless strain sensors, or fiber optic strain sensors. However, for this 
technical publication, surface strains were analytically calculated from finite-element analysis to 
demonstrate and validate the proposed structural deformation prediction methodology. 

By entering the surface strain data into the Displacement Transfer Functions, one can calculate slopes 
and deflections along each strain-sensing line on a given structure such as a wing. By using multiple 
strain-sensing lines, overall deformed shapes of a structure subjected to bending and torsion loadings can 
then be obtained. A total of seven sets of Displacement Transfer Functions were formulated in the past for 


different structural geometries (refs. 1-9), and were found to be quite accurate in the shape predictions of 
actual flight vehicles (refs. 10 and 11). 

By embodying the Displacement Transfer Functions, the rectilinearly distributed surface strains can 
also be input into the Stiffness and Load Transfer Functions to calculate structural stiffness (bending and 
torsion) and operational loads (bending moments, shear loads, and torques) for monitoring the operational 
loads of a flight-vehicle (ref. 12) in near real-time. The accuracy of this method for estimating operational 
loads on structures was analytically confirmed by using finite-element analysis of different aerospace 
structures such as tapered cantilever tubular beams, depth tapered un-swept wing boxes, depth tapered 
swept wing boxes, and the doubly-tapered Ikhana aircraft wing (ref. 13). 

All the earlier Displacement Transfer Functions (refs. 1-9) were formulated based on the Shifted 
straight deflections perpendicular to the undeformed neutral axis. For the shape predictions of structures 
under large geometric nonlinear deformations, there was a need to reformulate a new set of Displacement 
Transfer Functions for structure shape predictions. 

This technical publication is based on the actual curved displacement (curved distance traced by a 
material point from its undeformed position to its deformed position) and formulate the Curved 
Displacement Transfer Functions for shape predictions of structures under large geometric nonlinear 
deformations. 

A long tapered cantilever tubular beam was chosen to assess the accuracy of the new Curved 
Displacement Transfer Functions. Linear and nonlinear finite-element analyses were performed on the 
tapered cantilever tubular beam to 1) analytically generate linear and nonlinear surface strains for use in 
the Curved Displacement Transfer Functions and 2) analytically generate beam deflections (benchmarks) 
for comparisons with theoretically predicted beam deflections. The Curved Displacement Transfer 
Functions were then programmed using the finite-element-generated surface strains as input data to 
calculate theoretical deflections. The shape prediction accuracies were then determined by comparing the 
theoretical deflections with the finite-element-generated deflections (benchmarks). The results show that 
the Curved Displacement Transfer Functions, just like the Shifted Displacement Transfer Functions (refs 
1-9) are very accurate for shape predictions of structures under both linear and large geometric nonlinear 
deformations (beam-tip deflection reaching 58 percent of span). 


REVIEW OF TECHNICAL BACKGROUND 


To formulate the displacement theory (refs. 1-8), strain-sensing stations (strain measurement points) 
are to be discretely distributed along a strain-sensing line on the surface of the structure (for example, 
aircraft wing) (fig. 3). The depth-wise cross section of the structure along the strain-sensing line is called 
an imaginary “embedded beam” (not to be confused with the classical isolated Euler-Bernoulli beam). 
Each embedded beam was then discretized into multiple domains (strain-sensing station separation 
distances) with domain junctures matching the locations of the strain-sensing stations. By discretization, 
the variation of the embedded beam depth factor can be described with a piecewise linear function, and 
the surface strain variation can be described with either a piecewise linear or a piecewise nonlinear 
function. This approach enables the integration of the curvature equation of the deformed embedded beam 
to yield slope and deflection equations in recursive forms. The recursive slope and deflection equations 
are then combined into a single deflection equation in dual summation form. A set of three equations 
(recursive slope equation, recursive deflection equation, and dual-summation deflection equation) are 
called Displacement Transfer Functions, which are expressed in terms of the embedded beam geometrical 
parameters and surface strains, and contain no material properties. By entering surface strain data into the 
Displacement Transfer Functions, one can calculate slopes and deflections along the embedded beam. By 
using multiple strain-sensing lines, deflections at multiple strain-sensing stations can then be calculated 
for plotting the overall deformed shapes of the structure subjected to bending and torsion loads. The 
Displacement Transfer Functions are purely geometric in nature, and therefore, for a given density of 
strain-sensing stations, one can accurately compute the associated deflections, whether the input surface 
strains come from linear or nonlinear deformations. 


BASIC EQUATIONS FOR THE DISPLACEMENT THEORY 


The following first section is to geometrically establish the basic curvature-strain equation, from 
which different Displacement Transfer Functions can be formulated. The second section discusses 
different curvature-strain differential equations. 


Curvature-Strain Relationship 


Figure 4 shows the deformed state of an embedded beam with a changing depth factor, c(s), where s 


is the curved coordinate along the elastic curve of the deformed embedded beam. The curvature-strain 
relationship can be established graphically from figure 4. The embedded beam elastic curve has a local 
radius of curvature, R(s), within a small beam segment subtended by d@ . The undeformed curve length, 


AB, lies on the beam neutral axis, and the deformed curve length, A’B’ {= AB [1+ €(s)]} , where €(s) is 


the surface strain, lies on the beam lower surface. From the two similar slender sectors, O’AB and 
O’A’B’ , one obtains the geometrical relationship described with equation (1): 
O’A’ | a c(s) _ A’B’ 
O’A Ris) AB 


=1+e(s) (1) 


From equation (1), one obtains the curvature-strain equation (2): 


1 _ es) 
R(s) 7 c(s) 2 


Equation (2) geometrically relates the local curvature, 1/ R(s), to the associated surface strain, €(s), 
and the depth factor, c(s) , of the embedded beam. Equation (2) is the basis for formulating Displacement 
Transfer Functions. 


Traditional Curvature Equations for Vertical Deflections 


Different forms of curvature-strain differential equations written in the x-ysystem have the 


following familiar forms as shown in equations (3) through (5) (for detailed discussions, see refs. 14 and 
15). The mathematical differences between those equations are discussed. 


1, Eularian Curvature Equation 


1 2 da’ y/dx? _ €(x) 
R(x) [1+ (dy/dxY FP? e(x) 


(3) 


It is important to mention that equation (3) is referenced to the deformed (movable) x-coordinate (that is, 
x - coordinate gives only the deformed location of a material point, but not the undeformed location). 


2. Lagrangian Curvature Equation 


ie d°y/dx _ €(x) 


RX) fi—(dyjdxy (2) “ 


Equation (4) is in reference to the undeformed (fixed) x-coordinate (derivation in Appendix A, ref. 15). 


Because of the nonlinear term, (dy/ dx)’, direct integrations of equations (3) and (4) can end up in 


extremely complex deflection equations, which have poor prediction accuracies at large deformations 
(details in ref. 15). 


3. Shifted Curvature Equation 


If the deformed material points are shifted back to their respective undeformed x-positions [that is, by 
setting axial displacement u to zero (u— 0) ] (fig. 5), the shifting condition (u—0)_ will cause the 


nonlinear term, (dy/ any to become zero (Appendix A, ref. 15). Hence equation (4) becomes equation 


(5): 


fy EC) (5) 
dx” c(x) 

Equation (5) is in reference to the undeformed x-coordinate, is called the shifted curvature-strain 
differential equation, and is not the traditional linearized form of the Eulerian curvature equation (3) 
which is in reference to the deformed x-coordinate. 

Equation (5) is the basic curvature-strain differential equation used in the formulations of the Shifted 
Displacement Transfer Functions (ref. 15). Detailed discussions of equation (5) are provided in the 
subsequent Similarity of Shifted and Curved Formulations section. 


FORMULATION OF THE SHIFTED DISPLACEMENT TRANSFER 
FUNCTIONS 


The Shifted Displacement Transfer Functions were formulated earlier (refs. 1 and 4) by piecewise 
integrations of equation (5) for nonuniform embedded beams. The resulting equations are listed below as 
equations (6a), (6b), and (6c) and equations (7a), (7b), and (7c). As shown, equations (6a), (6b), and (6c) 
and equations (7a), (7b), and (7c) can degenerate into the limit cases of uniform embedded beams 


(c,, =c, =c). The limit cases were obtained by first expanding the logarithmic function, log,(c,/c,_,), 
in the neighborhood of (c,/c,,)=1, and then setting, c,, =c,=c. Equations (6a), (6b), and (6c) and 


equations (7a), (7b), and (7c) listed below were used in the linear shape prediction analysis portion of the 
present technical publication. 


Shifted Displacement Transfer Functions 


The following Shifted Displacement Transfer Functions shown as equations (6a), (6b), and (6c) were 
formulated by using the piecewise-linear representation of the variation of the depth factor, c(x), and 


using the piecewise-linear representation of the variation of the surface strain, €(x). The detailed 
mathematical derivations can be found in refs. 1 and 4. 


Slope equation: 


E, 7 EEC; — : 
tan@, = Ajj] "++ wie “log, “i fete, 
Ga-G (64-4 y C4 (6a) 
Uniform Al fa) 
—— 50 EH +e) +tane,, 
(i= 1,2,3,...,n) 
Vertical deflection equation: 
a. In recursive form: 
= €; é; c : 
= (Aly?| 27%: __ ert log, + +(c,,- «| +y,, +Altan@,, 
2(¢.4-C)) (C4 =a il 
(6b) 
Al 
— (Aly _, + Altan@,, 
Cj_j=C;=C 6c 
(i= 1,2,3,...,n) 
b. In summation form [equations (6a) and (6b) combined]: 
j E,4C;—E,C C; 
y, = (AL? > mst) ¢ log, + +(c),-¢,) 
ve —C, = (Cg ey Cr 
Contributions from deflection terms 
SE,  €3:56. -€.€6 
+ (Aly 3 ree fee a a 2 $= gg * —L |b+y,+@Altané, é 
ja Cy; (Cpe) Cj Sie se) 
=0 for cantilever beams 


Contributions from slope terms 
27 \ 
Uniform (Al ) 
SO ; f 
(¢j_)=C;=c) > (2€,_, + €j Ie 
6c “a 


Contributions from deflection terms Contributions from slope terms 


i-l 
(i— jE, + €;) +¥ + @Altané, 
1 


Je =0 for cantilever beams 


(i=1,2,3,...,2) 


Equations (6a), (6b), and (6c) are called the Shifted Displacement Transfer Functions for nonuniform 
embedded beams (c,,#c;), including the limit cases of uniform embedded beams (c,_, =c; =C) 


(refs. 1 and 4). 


Improved Shifted Displacement Transfer Functions 


The following Improved Displacement Transfer Functions shown as equations (7a), (7b), and (7c) 
were formulated for nonuniform embedded beams by using piecewise-linear representation of the 
variation of the depth factor, c(x), and using piecewise-nonlinear (quadratic) representation of the 
variation of the surface strain, €(x) (ref. 7). 


Slope equation: 


Al ic 
aie. = ae [ (2c, - G4 MCE — 26,4€)) + C,6,-4€ x1 Jlog, 2 : 

i i-1 i-l 

Al 
Fe [(Sc, —3¢,_,)€4 — 2c, —¢..4)€, + (c, +64)E 4, |+tan@_, (7a) 

i i-l 

Al 
= — (5€,_, + 8€; — €,,) + tan8,_, 
=O" 196 


(F122, 3s1058) 
Vertical Deflection Equations: 


a. Inrecursive form: 


(Al)? m 
= ————__| (2c. -c._, (c,€,_, —2¢,_,€.)+¢,.c._,€,,, || c, log, —-—-(c, -—c._ 
Yi 2(c,-¢,,)* [( i alt iVi-l i-l i) i~i-l val i Se 6s ( i iy) 
A 2 
~ AY 1 Be, = 5646-4 —2(5c, —2e,.)8, + (26, FC JEuy]+¥.q + Altand,, (7b) 
12(c,-—¢,4) 
Al)’ 
— (A) (7€,_, + 6€, —€,,,)+y,_, + Altan@,, 
eae 24c 


(i= 1,2,3,...,n) 
b. In summation form [equations (7a) and (7b) combined]: 


1 


i |%c.-c..)4 
y,= (any oH? 
j=l 


[@e, —€,4NC£,)—2€,,€,) +c 


iG: 
J 
i€1€ jan || C; log. : —(€;-¢,4) 


j-l 


a [ (8c, -— Se,_,)€,4—-2(Se, —2c, JE, + (Ze, Fe, )E pr | 
i j-l 


Contribution from deflection terms 


1 c, 
= ye Cp MEE py 2Op EFCC 48 ja Jog. Cc, 
+(AlIY SH G- jy x 
- (70) 
j arrears = 362 )64 = 2136 '6,5 EFC, +0; Ej | 


Contributions from slope terms 
+ y, +i(Al)tan8, 
ee 


=0 for cantilever beams 


Uniform (Al)° (Al) \ . . . 
ae yr pt 6€;—En,)+ cre — j(S5€,,+8€,-€),)+ yy + @Altan@, 
im iA =0 for cantilever beams 
Contributions from deflection terms Contributions from slope terms 


(i =1,2,3,...,n) 


Equations (7a), (7b), and (7c) are called the Improved Shifted Displacement Transfer Functions for 
nonuniform embedded beams (c,,#c,), including the limit cases of uniform embedded beams 


(6 Sc) Get 7), 


FORMULATION OF CURVED DISPLACEMENT TRANSFER 
FUNCTIONS 


The following sections present mathematical processes needed for the formulation of the new Curved 
Displacement Transfer Functions. The formulation is based on curved deflections instead of traditional 
vertical deflections. 


Curved Curvature Equations 


For large bending deformations of beams (fig. 5), one must understand that the actual (true) 
deflection, y, of a material point is a curved distance traced by the same material point from its initial 


undeformed position to its final deformed position. Thus, the conventional deflection, y, is merely the 
vertical component of the curved deflection, y . Also see discussions about large deformations in refs. 16 


and 17. 
The basic curvature equation referenced to the curvilinear s - y system, instead of traditional 


Cartesian x - y system, can be expressed as equation (8): 


1 do d(dy\ dy dy 
0). 4(8). 5, ga ‘ 


R(s) ds ds\ ds) ds? 


Equating equations (2) and (8) gives the curvature-strain differential equation (9) in the s- y system 
for large deformations: 


d*y _ &(s) 
ds? c(s) 


(9) 


Equation (9) is a purely geometrical relationship, containing no material properties. Assuming the 
length of neutral axis of the embedded beam remains the same (that is, 5s = x ) after bending, equation (9) 
can be rewritten in reference to the undeformed x-system as equation (10): 


dx’ c(x) (10) 


The mathematical process for formulating the Curved Displacement Transfer Function is through the 
piecewise integration of equation (10) and is described as follows. 


Piecewise Representations 


To enable piecewise integrations of equation (10), the depth factor, c(x), and the surface strain, 
€(x), can be expressed by either piecewise linear or piecewise nonlinear functions, described as follows 
(fig. 3) (refs. 1-8). 


1. Depth factors 


The variations of the embedded beam depth factor, c(x), within each small domain, x,_, <x <x; 
(i =1,2,3,...,2) , can be expressed with the linear function given by equation (11): 
T= 55 
(4) = 614 +(€, 6.) 4 : (%,S2=%,) (11) 


2. Surface Strains 


The variation of the surface bending strain, €(X), within each small domain, X,, [xX <xX;, can be 


expressed by either a linear function given by equation (12) (ref. 4) or by a nonlinear function given by 
equations (13a) and (13b) (ref. 7): 


a. Linear: 


EQ) = 6.7 +(,-E.. b 434 SaHSe.) (12) 
b. Nonlinear: 
3e€.,-4é€.+€ QE. +€ 
e(x)=e,, -— xx + Set 28 Fiat Sx SX; 
(x)= E, ay (4— 24) TINT (x-x,,)° ew ‘) (13a) 
6 S64 = 36,41 3, : (at i=n) (13b) 


Equation (13a) was generated by standard quadratic interpolation of strain values, {€, ,,€,,€,,,}. 
respectively at three equally spaced strain-sensing stations, {x,_,,x,,x,,,}, and equation (13b) is the 
quadratic extrapolation equation to obtain extrapolated strain, €,,, , beyond the embedded beam tip (p. 33, 


ref. 7). 


i+1? 


Piecewise Integrations 


In view of equations (11), (13a), and (13b); the curvature-strain differential equation (10) can be 
piecewise integrated to yield the Curved Displacement Transfer Functions. The piecewise integration of 
equation (10) within the domain, x,, < x <x,, between the two adjacent strain-sensing stations, {x,_,,x,} 


i-| — 


, yields the slope-angle equation (14): 


. Ps . * : 
J Trax 2/2) =a a5, ‘ Gi5 S425,) 
Xi-1 ax, ax, dx i-l Xj-1 c(x) (14) 
Eq. 6) a(x) 4x Eq.) 


i-l 


which can be rewritten in the form of equation (15): 


SD 


0 —dx+ 0. > (x;_ SxS X;) 
@=]o 6,4 AG 
slope 
Slope increment at x; 
Integration of the slope angle equation (14) yields the curved deflection equation (16): 
x | dy (dy 2 - EX) 
i ae WN ape x=y(x)- Yy ots G.de= J. ai 
x) dx \dx J. oa, te dicey (16) 
A(x) = at X;_| Deena at x ‘Right hand side. hand gide 
- due to 6, of equation (12) 
is == x) 
In view of equation (15), equation (16) can be rewritten in the form of equation (17): 
2 E(x) 
y(xX)= [ [ ——_ Wide PX Oa Ft Vex =f. Oa)det ¥, i 
Xj 9 Xj4 C(x) aE oa Kw x 
Deflection at x Deflection (17) 
Deflection increment above y;_; due to 6; at x} 


(Ag = 2S I.) 


t 

Using piecewise representations of {c(x),€(x)} given by equations (11), (13a), and (13b); equations 
(15) and (17) can be integrated within the domain, x,,<x<4xX;, to yield the slope and deflection 
equations in closed recursive and summation forms. A set of three equations (recursive slope equation, 
recursive deflection equation, and summation deflection equation) are called the Curved Displacement 
Transfer Functions. The mathematical processes are similar to those used in the piecewise integrations of 
the shifted curvature-strain gifferential equation (5) to formulate the Shifted Displacement Transfer 
Functions (refs. 1, 4, and 7). The detailed mathematical derivations of the Curved Displacement Transfer 


Functions are presented in Appendices A—H. 


CURVED DISPLACEMENT TRANSFER FUNCTIONS 


After piecewise integrations of equations (15) and (17), and going through a mathematical process 
ainiiag © the one used in formulating the Shifted Displacement Transfer Functions (refs. 1, 4, and 7), one 
can obtain the Curved Displacement Transfer Functions shown as equations (18a), (18b), (18c), (19a), 
(19b), and (19c) for nonuniform embedded beams, including the limit cases of uniform embedded beams, 


(c_, =¢,=Cc). 


Curved Displacement Transfer Functions 


The Curved Displacement Transfer Functions shown in equations (18a), (18b), and (18c) were 
formulated by carrying out piecewise integrations of equations (13) and (15) using the piecewise-linear 
representation of the depth factor c(x) [eq. (11)] variation, and using the piecewise-linear representation 


of the surface strain €(x) [eq. (12)] variation. The detailed mathematical derivations are presented in 
Appendices A, B, C, and D. 
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Slope-angle equation (Appendices A and C): 


Ero E:. E : 
0, = Al i-l i ap i- -1Cj fel log, C; at 60. , 
C1 7 Gj (¢4- C;) Ci 
niform Al (18a) 
— 30 eH +€)+0,, 
(i= 1,2,3,...,n) 
Curved deflection equation: 
a. In recursive form (Appendices A and C): 
—€, E. Eg, ; a 
(A) a cae ee Ce 7) +5, +(ADO,, 
AGC) (C4 —¢,)° il 
(A di (18b) 
Pee eu te) Aes 
(i= 1,2,3,...,n) 
b. In summation form [equations (18a) and (18b) combined] (Appendices B and D): 
j Ex Cj E:C Cc. 
y, =(Al) Dees eee ae A Clg. F(t) 
Te =e . (c, -c,) G4 
Contributions from deflection terms 
le €,-€ €jC.-ECc ic 
+ (Al) ¥°) (i- p) fate —E i gg “| + Jy + (ANB, 
Jel Cj 7 © (Cj ~ cy Cj =0 for cantilever beams qd 8c) 


Contribution from slope terms 
2 i-l 


vito cay GY 0c, 6) +S Ql DE +E,) + Fo + O(ADB, 


(c,4=¢;=¢) 


=0 for cantilever beams 


Contributions from deflection terms Contributions from slope terms 


(i=1,2,3,...,n) 


Equations (18a), (18b), and (18c) are called the Curved Displacement Transfer Functions for nonuniform 
embedded beams (c,#c,,) (Appendices A and B) under large deformation with geometrical 


nonlinearity including the limit cases of uniform embedded beams (c, , = c, =c) (Appendices C and D). 


Improved Curved Displacement Transfer Functions 


The Improved Curved Improved Displacement Transfer Functions shown in equations (19a), (19b), 
and (19c) were formulated by carrying out piecewise integrations of equations (13) and (15) using a 
piecewise-linear representation of the depth factor c(x) [eq. (9)] variation, and using a piecewise- 
nonlinear representation of the surface strain e€(x) [eqs. (1la) and (11b)] variation. The detail 
mathematical derivations are presented in Appendices E, FP, G, and H. 
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Slope-angle equation (Appendices E and G): 


Al C, 
0. = ———_| (2c. -c._, (c.€._, —2c,_,€,)+c.c,_,€.,, log, — 
i 2(c, e ci) [( i i-l i~i-l i-l i) i-i-l cal Se Ch 
Al 
cae [(Sc;— 3¢,.1)€)-1 — 2(3¢; — 8; + (C+ 6-1)E ina |+ 8 (19a) 


; Al 
( 7 — (964 a 8E, = Ei) 6. 
Mark De 


(i=1,2,3,....n) 


Curved-deflection equations: 


a. In recursive form (Appendices E and G): 


2 Al)’ . 
= _ AN Ie, — C4 MCE — 26 4€)) + CC Ei I< log, . =(O C5 ) 


= 4 
—€,4) 74 


yi 2c, 


Aly - 
a Delite — 5c), )E — 2(5c, — 2¢,,)€; + (2c, +64 Ein |+. 54+ AO, (19D) 
i i-l 

Uniform (Al y 


es ay A + 6€,-€,,)+y,_, +(AD@,, 


(i=1,2,3,...,n) 


b. In summation form [equations (19a) and (19b) combined]: (Appendices F and H): 


1 C. 
i yee le — Cj MCE p44 —2€; E+ G0) 52 i | clog, ———(c, C4) 
= 2 (c;-¢;4) Cp 
aC» 
j=l 
“Tie me gL es Sees — 21565 — 26, Hey HEE 


Contributions from deflection terms 


1 cj 
7 ——— > i@e, C24 MCE 2698) FCC 8 2, log. — 
Haney Gp) a 
j=l 
(Se; — 3c; )€;1—2Be;—¢, Je; +(e, FEE ja | 


ee 


Contributions from slope terms 
+ yy + ()(Al)O, 
=0 for cantilever beams 
niform (Al) (Al) — : . = : 
eas ve YS (Te, +6€, Ej) + oe YG = PSE + 8E; — Ej) + Yo + DADO, 
jel = 


Contributions from deflection terms 


(19c) 


=0 for cantilever beams 


Contributions from slope terms 


(i=1,2,3,....n) 
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Equations (19a), (19b), and (19c) are called the Improved Curved Displacement Transfer Functions 
for nonuniform embedded beams (c;#c,,) (Appendices E and F), under large deformations with 


geometrical nonlinearity including the limit cases of uniform embedded beams (c,, =c, =c) 


(Appendices G and H). Equations (18a), (18b), (18c), (19a), (19b), and (19c) listed above were used in 
the nonlinear shape prediction analysis portion of the present technical publication. 
It is important to mention that if {y,,0,} in equations (18a), (18b), (18c), (19a), (19b), and (19c) are 


replaced respectively with {y,,tan0,}, then the Curved Displacement Transfer Functions shown in 


equations (18a), (18b), (18c), (19a), (19b), and (19c) will become the Shifted Displacement Transfer 
Functions shown in equations (16a), (16b), (16c), (17a), (17b), and (17c). As will be seen shortly, the 
vertical (straight) and curved deflections, {y,, y,}, calculated respectively from the Shifted and Curved 


Displacement Transfer Functions turned out to be identical. 


CHARACTERISTICS OF DISPLACEMENT TRANSFER FUNCTIONS 


In the Shifted and Curved Displacement Transfer Functions, the vertical and curved deflections, 
{y,,¥;} , at the strain-sensing station, x,, are expressed in terms of the inboard depth factors 
(Cy .C),C>,+.5C;) and the associated inboard surface strains (E, ,€,,€,-...€;) [for eqs. (6a), (6b) (6,c); eqs. 
(18a), (18b) (18c)] or (€,€,,€5,...5€;,,) [for eqs. (7a), (7b), (7c); eqs. (19a), (19b), (19c)], including the 
values of {c,,€,} at the strain-sensing station, x,, where deflections, {y,,,;}, are calculated. 

It is important to mention that equations (6a), (6b), (6c), (7a), (7b), (7c), (18a), (18b), (18c), (19a), 


(19b), and (19c) are purely geometrical relationships, containing no material properties. However, it must 
be understood that the surface strains, €,, implicitly contain the effect of material properties and internal 


structural configurations. Thus, in using equations (6a), (6b), (6c), (7a), (7b), (7c), (18a), (18b), (18c), 
(19a), (19b), and (19c) for shape predictions of complex structures such as aircraft wings, there is no need 
to know the material properties, nor the complex geometries of the internal structures. 


STRUCTURE USED FOR SHAPE PREDICTION ANALYSES 


The structure chosen for shape prediction analysis is a long tapered cantilever tubular beam with 
dimensions and material properties listed in table 1. 


Table 1. Dimensions and material properties of a tapered cantilever tubular beam. 


Beam root dimensions Beam tip dimensions Material properties 
(for Nastran analysis) 


1, in t, in hy in Co» in h Be 0) E, lb/in* Vv 


300 0.02296 8 4 ) 1 10.5 x 10° 0.3 


Figure 3 shows a span-wise vertical cross-section (embedded beam) of the tapered cantilever tubular 
beam with strain-sensing lines located on both lower and upper surfaces. The embedded beam is 
discretized into n = 20 equal domains with domain length of Al =1/n= (300 in.)/20 =15 in. Thus, there 
are 1+1=21 equally spaced strain-sensing stations along each strain-sensing line. Different lateral loads P 
={50, 100, 200, 300, 400, 500, 600}Ib were applied at the beam tip to study the effect of geometric 
nonlinearity on the shape predictions under increasing loads (deformations). For the tubular beam, the 
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local depth factors, c, ((=1,2,3,...,n), are the local outer radius of the tubular beam, and are known. As 
will be shown in the finite-element analyses, for the linear cases, the depth factors, c,, remain unchanged 


because, the magnitudes of the pairs of lower and upper surface strains are identical regardless of the load 
level. Therefore, only the lower surface strains are needed for inputs to the Displacement Transfer 
Functions for shape calculations. 

However, for large geometric nonlinear deformations (for example, the Helios flying wing shown in 
figure 1), the lower and the upper strain-sensing lines are needed because the magnitudes between the 
lower and upper surface strains will be slightly different due to the axial strain components induced by the 
curved-beam effect at large bending. Thus, both lower and upper surface strains are needed to calculate 
the deformation-dependent depth factors, c,, which together with associated lower surface strains, €,, 


can then be input to the Displacement Transfer Functions for shape calculations. 


FINITE-ELEMENT ANALYSES 


The MSC/Nastran (MSC Software Corporation, Newport Beach, California) finite-element program 
(ref. 18) was used for both linear and nonlinear analyses of the tapered cantilever tubular beam subjected 
to different beam-tip loads. In this technical publication, the surface strains needed for input to the 
Displacement Transfer Functions for shape predictions were Nastran-generated and not experimentally 
measured. Also, the Nastran-generated deflections were used as reference benchmarks to study the shape 
prediction accuracies of the Displacement Transfer Functions. 

Figures 6a and 6b show two types of finite-element models generated for the tapered cantilever 
tubular beam using shell elements (fig. 6a) and using beam elements (fig. 6b). When modeling the tubular 
beam with the shell elements, the nonlinear analysis could be carried out only up to certain low level 
nonlinear deformations, beyond which the nonlinear analysis would break down. Namely, the shell 
element stiffness would become ill conditioned (that is, mathematical singularity) due to a very large ratio 
between bending and membrane stiffness in a large displacement situation. When the beam elements were 
used, the nonlinear analysis could be carried out up to very large deformations without encountering the 
above-mentioned mathematical singularity. 


Nastran Linear Analysis 


The Nastran linear analysis (using a linear strain tensor) assumes a linear relationship between the 
load applied to a structure and the response of the structure. In using linear theory for large deformation 
analysis, the deflection of the structure is simply proportional to the apply load. 

In the linear analysis, the Nastran displacement outputs provide vertical deflections, y , but zero axial 
x-displacements (u=0) (fig. 6a). The zero axial displacement implies the horizontal (x-direction) 
shifting of the vertical deflections to their respective undeformed x -positions. Therefore, the vertical 
deflections generated by the Nastran linear analysis can be compared with the vertical deflections 
calculated from the Shifted Displacement Transfer Functions using the Nastran-generated linear surface 
strains. 


Nastran Nonlinear Analysis 


In the geometric nonlinear large deflection problem (strain tensors contain second order terms) the 
stiffness of the structure depends on the displacement, and thus the deflection response is no longer a 
linear function of the load applied (refs. 18 and 19). For nonlinear analysis, both x- and y- 


displacements, {u,y}, are nonzero along the neutral axis (fig. 6b and fig.5). In the large deformations, 


the load is no longer vertical; it has followed the structure to its deformed state (fig. 6b). In the Nastran 
nonlinear analysis, the Follower Force Option command was used to force the applied load to be normal 
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to the deformed neutral axis without introducing an axial load component (fig.6b). For comparisons with 
the curved deflections, y, predicted from the Curved Displacement Transfer functions, the Nastran- 


displacement data of {u,y} were used to generate the equivalent curved deflections, y, by using the 
deflection-conversion equations described in the following section. 


STRAIGHT-TO-CURVED DEFLECTION CONVERSION 


As previously mentioned, the Nastran linear analysis provided only vertical displacements, y, but 
zero axial displacements, (u =0), along the neutral axis (fig. 6a). However, in the Nastran nonlinear 
analysis, the Nastran outputs gave both vertical displacements, y, and axial displacements, u(#0), 
along the neutral axis (fig. 6b). To compare with the curved deflections, y, calculated from the Curved 
Displacement Transfer functions, the Nastran nonlinear displacement data of {u,y} were used to 
generate Nastran versions of curved deflections, y. 

Figure 7 shows the elastic curve of the deformed embedded beam. In figure 7, 
6{= tan'[y,/(x, —u,)]}, is the slope angle of a slanted straight line connecting the origin and deformed 
point, s;, on the embedded beam elastic curve, and u,; is the magnitude of the axial displacement of 
deformed material point, s,, along the undeformed x -axis. In figure 7, if the curved deflection, y, is 


considered as a circular arc, then 6. will be an angle subtended by the arc length, y, which can then be 


related to the vertical deflection, y, , through the following deflection-conversion equation (20) (fig. 7). 


y= x0, =X; tan” * 
Lh (20) 


Equation (20) was used to convert Nastran displacement data of {u,,y,} into Nastran versions of curved 


deflections, y,, for comparison with theoretically predicted curved deflections, y,, from the Curved 
Displacement Transfer functions. 


CURVED-BEAM EFFECT IN NONLINEAR DEFORMATIONS 


In the Nastran linear finite-element analysis of the tapered cantilever tubular beam, the Nastran strain 
outputs showed that the magnitudes of the lower and upper bending strains at the same cross section are 
always identical regardless of loading levels. The equal magnitudes of the lower and upper bending 
strains implies that the neutral axis is always located at the half depth of the embedded beam, and no axial 
strains are induced under linear bending. 

However, in the Nastran nonlinear analysis, the Nastran strain outputs show that the magnitudes of 
the lower and upper surface strains at the same cross section are not exactly the same under large bending. 
The magnitudes of the lower surface strains are slightly larger than the magnitudes of the associated upper 
surface strains especially in the outboard flexible region. The magnitude difference of the lower and upper 
surface strains could be attributed to a slight shifting of the neutral axis toward the bent concave side of 
the embedded beam due to curved-beam effect. 

Before using the Displacement Transfer Functions, the correct depth factor and surface strains must 
be used. The following two methods: 1) Neutral-axis shifting method or 2) Axial-strain elimination 
method can be used to calculate the correct depth factor or correct surface bending strains for use in the 
Displacement Transfer Functions. 
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Neutral-Axis-Shifting Method 
As shown in figure 8, if the difference between the lower and upper surface strains magnitudes is due 
to neutral axis shifting, then one can use pairs of the lower and upper strains, {€,,€,} , to calculate the 
unknown lower depth factors, c; (or upper depth factor c,)((=1,2,3,...,.n), from the depth-factor 
equation (21): 


C= ih, : C.=h.-c. : (i =1,2,3,...,n) (21) 


in which h, (i =1,2,3,...,7) is the depth of the embedded beam at x = x,. For the current linearly tapered 
embedded beam (fig. 3), the depth, 4, , can be calculated from the following depth equation (22): 


h, = Ny — “(hy “hy = 21793000) (22) 


in which {h,,h,}are respectively the depths of the embedded beam at the embedded beam root 
(x = xX, =O) and at the embedded beam tip (x =x, =/). 


When the shifted lower depth factor, C; (eq. 21), is used, the associated lower surface strain, €, , must 


also be used for input to the Displacement Transfer Functions. As will be seen shortly, the neutral axis 
shifting method can automatically nullify the axial strain effect. 


Axial-Strain Elimination Method 


If the known depth factor c, =h,/2 of the embedded beam (fig. 3) is to be used, one can consider the 


unequal magnitudes of the lower and upper surface strains contain both bending and axial strain 
components. In view of figure 9, by averaging the magnitudes of the lower and upper surface strains 
(€, >0,€, <0), axial strain components can be eliminated to yield the true bending strains given by 


equation (23) (fig. 9): 


, ; €,—€, 
True bending strain = =e ‘ (i = 1,2,3,...,.7) (23) 


For entering the true bending strains, (€, — €,)/2 , into the Displacement Transfer Functions, the known 
depth factor, c,(= h, /2), must be used. 


Strain-to-Depth Factor Ratios 


For exploration purpose, the strain-to-depth factor ratios for the above two methods [eqs. (21) and 
(23)] can be written in the following forms as equations (24) and (25): 


For neutral-axis-shifting case [rewriting eq. (21)]: 


(24) 
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For axial-strain elimination case [from eq. (23)]: 


H( £58.) -2( £58) 2-8 25 
cl 2 ) al 2 h, -) 


Note from equations (24) and (25) that the strain-to-depth factor ratios for the two methods are 
identical. Since all the Displacement Transfer Functions from equations (6a), (6b), (6c), (7a), (7b), (7c), 
(18a), (18b), (18c), (19a), (19b), and (19c) are expressed in terms of strain-to-depth factor ratios, the 
slopes and deflections calculated using the two input methods turned out to be extremely close as will be 
seen in the Numerical Results section. 


PREDICTION ERROR EQUATIONS 


The Nastran-generated deflections were used as benchmark data to study the theoretical deflection 
prediction errors. Let {y,,y,}, respectively denote the vertical and curved deflections predicted, 


respectively from the Shifted and the Curved Displacement Transfer Functions, and let {(y,)y,.(¥)yw}> 


respectively denote the corresponding Nastran-calculated linear and nonlinear deflections. Then, the 
prediction error is defined by the following prediction error equations (26) and (27), respectively for 
vertical and curved deflection cases: 


Vertical deflection case (linear analysis): 


Prediction error = [2 - J x 100% (26) 
(Y) ne 


Curved deflection case (nonlinear analysis): 


Ji) NN 


Prediction error = Me ] x 100% (27) 


Equations (26) and (27) were used to determine the prediction errors of the Shifted and Curved 
Displacement Transfer Functions. 


NUMERICAL RESULTS 


The complete set of strain and deflection data generated by Nastran linear and nonlinear analyses of 
the tapered cantilever tubular beam are tabulated in Appendix I for different beam-tip load P. 


Nastran Linear Cases 


All the data generated for the Nastran linear cases are listed in tables Il1-I7 of Appendix I. Note from 
tables I1—I7 that for the Nastran linear cases, the lower and the upper surface strains at the same strain- 
sensing cross sections have the same magnitudes. The Nastran outputs gave only vertical deflections, y, , 
and zero axial displacements (that is, vu, =0) (fig. 6a). Thus, Nastran linear case is equivalent to the 
Shifted formulation. In the last columns of tables I1—I7, the corresponding theoretical vertical deflections, 
y,; , were calculated from the Shifted Displacement Transfer Functions in equation (6b) using the known 
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depth factors, c,(=h, / 2), and Nastran lower surface strains, €,, listed in tables I1-I7. The theoretical 


deflections, y,, calculated from equation (7b) are not listed because both equations (6b) and (7b) gave 
practically identical vertical deflections with maximum difference of only 0.08 percent at the beam-tip for 
the present low-tapered cantilever tubular beam. The data listed in tables I1-I7 of Appendix I were used 
in plotting the following linear strain curves and vertical deflection curves for visual display. 


I. Nastran Linear Strain Curves 


Figure 10 shows surface strain curves associated at different load levels generated from Nastran linear 
analysis of the tapered cantilever tubular beam based on data listed in tables I1— I7 of Appendix I. Note 
that for linear cases, the magnitudes of the lower and upper surface strains at the same axial location are 
identical, and increase linearly with increasing load, P. Note from figure 10 that for the present tapered 
tubular beam, with tip-to-root depth ratio, (c, /C)) = 1/4, the linear strains increase almost linearly in the 
span-wise direction in the inboard regions, reaching the peaks in the outboard regions, and then decrease 
rapidly down to zero at the beam tip. 


2. Vertical Deflection Curves 


Figure 11 shows vertical deflection curves for the tapered cantilever tubular beam at different loading 
levels calculated from the Shifted Displacement Transfer Functions [eq. (6b)] and from Nastran linear 
analysis based on the data listed in tables I1—-I7 of Appendix I. 

As shown in figure 11, the theoretical deflection curves [eq. (6b)] practically fell on top of the 
corresponding Nastran-generated deflection curves, even up to very large bending under P = 600 Ib, with 
beam-tip deflection reaching 94 percent of the beam span (see table I7 of Appendix I), and beam-tip slope 
angle reaching 66 deg. (ref. 15). The good agreement between the vertical (straight) deflections calculated 
from the Shifted Displacement Transfer Functions and from Nastran linear analysis, gives confidence in 
the mathematical formulations of the Shifted Displacement Transfer Functions. 


Nastran Nonlinear Cases 


All the data generated for the Nastran nonlinear cases are listed in tables [I8—I14 of Appendix I. The 
Nastran outputs gave both axial and vertical displacement, {u, ,y,} (fig. 6b). Thus, Nastran nonlinear 
analysis is equivalent to the Curved formulation. Note from tables I8—I14 that the magnitudes of the lower 
surface strains are slightly larger than the magnitudes of the associated upper surface strains. Such lower 
and upper strain magnitude differentials could be attributed to the curved-beam effect, which induces 
slight axial tensile strain components under nonlinear bending. 

In the last columns of, the theoretical curved deflection, y, were calculated from the Curved 
Displacement Transfer Functions [eq. (18b)] using the known depth factors, c,(=h,/2), and the true 
bending strains (€,—€,)/2 [eq. (23)] using the Nastran-nonlinear strain data, {€, ,€,},, listed in tables 
18-114 of Appendix I. The theoretical deflections, y,, calculated from equation (19b) are not listed 
because both equations (18b) and (19b) gave practically identical curved deflections for the present low- 
tapered cantilever tubular beam case. For highly tapered flexible beam cases, equations (18b) and (19b) 
can give slightly different curved deflections, y, (ref. 7). The data listed in tables I8—I14 of Appendix I 
were then used in plotting the following nonlinear strain curves and curved deflection curves for visual 
display. 
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1. Nastran Nonlinear Strain Curves 


Figure 12 shows surface strain curves associated with different load levels generated from Nastran 
nonlinear analysis of the tapered cantilever tubular beam using data listed in tables I7-I14 of Appendix I. 
In the plots of the upper surface strains, which are negative, only the magnitudes were used (dashed 
curves). 

For the nonlinear cases, the magnitudes of lower and upper surface strains at each load level are no 
longer identical. Note from figure 12, that for the load less than P=100 lb, the lower and upper strain 
curves are equal. However, when the load level exceeds P=100 Ib, the magnitudes of upper strains in the 
outboard region become slightly less than the corresponding lower strains, and the difference between 
each set of lower and upper strain curves increases with increasing load. The difference between the 
lower and upper strain magnitude for the present tubular beam case can be attributed to the curvature 
effect induced by nonlinear bending, causing slight neutral axis shifting toward the compression boundary 
(fig. 8). 


2. Curved Deflection Curves 


Figure 13 shows the curved deflection curves for the tapered cantilever tubular beam at different 
loading levels calculated from the Curved Displacement Transfer Functions [eq. (18b)] and from Nastran 
nonlinear analysis. For convenience, the horizontal displacements are neglected, and the deformed 
material points were plotted on their respective undeformed x-locations. 

As shown in figure 13, the theoretical deflection curves [eq. (18b)] practically fell on top of the 
corresponding Nastran-generated deflection curves, even up to very large nonlinear bending under P = 
600 Ib, with beam-tip deflection reaching 58 percent of the beam span, and beam-tip slope angle reaching 
69 deg. (see table 2 and table 114 of Appendix I). 

Lastly, by using Nastran nonlinear strains as inputs, the theoretical curved deflection curves 
[calculated from eq. (18b)] and the corresponding Nastran-generated curved deflection curves [calculated 
from eq. (20)] are graphically coincidental (fig. 13). The excellent agreement between the curved 
deflections calculated from the Curved Displacement Transfer Functions and from Nastran nonlinear 
analysis, gives confidence in the mathematical formulations of the Curved Displacement Transfer 
Functions. 


Comparisons of Neutral-Axis Shifting Method and Axial-Strain Elimination Method 


Based on Nastran-generated nonlinear surface strains for a typical case of P=600 Ib listed in table 114 
of Appendix I, both the neutral-axis-shifting method and the axial-strain-elimination method were used to 
obtain correct data for input to the Curved Displacement Transfer Functions [eq. (18b)] for calculations of 
the slope angle, @,, and curved deflections, y,, for nonlinear deformations. The resulting data are 


compared in table 2. 
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Table 2. Comparisons of slopes and curved deflections, {0,,y,} , calculated from Curved Displacement 
Transfer Functions [eq. (18a) and (18b) for the P=600 lb nonlinear case using two input methods: 
1) Neutral-axis-shifting method and 2) Axial-strain-elimination method based on the Nastran nonlinear 
strain data of table 114 of Appendix I. 


Neutral-axis-shifting method Axial-strain-elimination method 
Using calculated C, [eq.(21)] and Using known ¢,(=h,/2) and 
Nastran nonlinear lower surface strains E, (table I14) true bending strains (€; — E;) ij 2 [eq. (23)] 
i c;, if 6, , deg y,, in ¢,, in 6, , deg y,, in 
Calculated eq. (21) Eq. (18a) Eq. (18b) Known (= h,/2) Eq. (18a) Eq. (18b) 
0 4.002 0.000 0.000 4.000 0.000 0.000 
1 3.852 2.541 0.330 3.850 2.541 0.330 
2 3.703 5.208 1.342 3.700 5.208 1.342 
3 3.554 8.006 3.069 3.550 8.006 3.069 
4 3.406 10.940 5.546 3.400 10.940 5.545 
5 3.258 14.018 8.809 3.250 14.018 8.809 
6 3.109 17.243 12.898 3.100 17.243 12.898 
7 2.961 20.619 17.851 2.950 20.619 17.851 
8 2.814 24.149 23.708 2.800 24.149 23.708 
9 2.666 27.836 30.509 2.650 27.835 30.509 
10 2.519 31.675 38.296 2.500 31.675 38.296 
11 2.371 35.662 47.107 2.350 35.662 47.107 
12 2.224 39.786 56.981 2.200 39.785 56.980 
13 2.077 44.027 67.949 2.050 44.027 67.949 
14 1.931 48.355 80.041 1.900 48.355 80.040 
15 1.785 52.724 93.211 1.750 52.724 93.271 
16 1.641 57.058 107.644 1.600 57.057 107.643 
17 1.501 61.227 123.133 1.450 61.227 123.132 
18 1.367 65.013 139.669 1.300 65.015 139.668 
19 1.264 67.992 157.103 1.150 68.001 157.104 
20 _|_No data (0/0) 69.4394 175.138* 1.000 69.416" 175.138* 
ANegligible differences (0.0331%) * Identical at beam tip, 


Notice from table 2 that the slope angles, 0;, and curved deflections, y, , calculated respectively from 
equations (18a) and (18b) based on the Neutral-axis-shifting method and the Axial-strain-elimination 
method are extremely close. These results show that either method could be used for the present tubular 
beam case, for which the neutral axis is located at the half depth of the embedded beam. Remember that 
the Axial-strain-elimination method is used to eliminate axial strains only when the neutral axis is located 
at the half depth of the undeformed embedded beam. For a complex structure (for example, aircraft 
wings) with unknown neutral axis location, the Neutral-axis-shifting method can be used to calculate the 
unknown depth factors, c,. 

Figure 14 shows the plots of both calculated and known depth factors, c;, using the c,-data listed in 
table 2. Note that the calculated depth-factor curve lies slightly above the known depth-factor curve, 
c,(=h,/2), implying slight shifting of the neutral axis toward the concave (compression) side of the 
deformed embedded beam, especially in the outboard region. The reason for the neutral axis shift is that 
under large deformations, the difference between the magnitudes of the lower and upper surface strains at 
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the same cross section increases in the highly bent outboard region (fig.12). Note that at the beam tip, no 
data point is shown (table 2) because the beam-tip surface strains are zero, causing equation (21) to give 
indefinite 0/0 value. 


Linear-Nonlinear Transition 


Figure 15 shows the plots of Nastran-generated beam-tip deflections, {y,,y,}, as functions of 


applied load, P. The corresponding theoretical deflection curves calculated from the Shifted deflection 
equation (6a) and Curved deflection equation (18a), respectively using Nastran linear and nonlinear strain 
data, graphically fell on top of the corresponding Nastran deflection curves. Using Nastran linear strains, 
the vertical deflection, y,, is a linear function of applied load, P. Therefore, the linear large deformation 
is simply the scaled up version of the small deformation. However, using Nastran nonlinear strains, the 
curved deflection, y,, increases convex upwardly with the applied load, P. Note from figure 15 that the 


linear and nonlinear beam-tip deflection curves are practically the same up to P=100 lb, at which the 
beam-tip deflection to span ratio is (y, /1) =(y,,/1) = 0.156, Beyond P=100 Ib, the linear and nonlinear 


beam-tip deflection curves diverge. Therefore, for the present tapered cantilever tubular beam, the 
normalized deflection of (y, //) = 0.156 can be considered as the borderline between linear and nonlinear 


deformation regimes. 
For the current long tapered cantilever tubular beam under the tip load of P=600 Ib., the beam-tip- 
deflection-to-span ratio is (y, /l) = (y, /J) = (175.138 in.)/(300 in.) = 0.58, which is larger than the Helios 


case of (beam-tip-deflection)/(half-span) = (40 ft)/(123.5 ft) = 0.32. 


Similarity of Shifted and Curved Formulations 


Note that the curvature-strain differential equations (5) and (8) have identical right-hand sides, 
€(x)/c(x). For exploratory purposes, the vertical deflection, y, calculated from equation (6b) (Shifted 


formulation) and the curved deflection, y, calculated from equation (18b) (Curved formulation) are 
compared for the P=600 lb nonlinear case. In the calculations of deflections, {y,y}, the calculated depth 
factors, c; [eq. (21)], and the nonlinear lower surface strains, €, (table 114 of Appendix I), were used as 
inputs. Table 3 lists the calculated results for the P=600 Ib nonlinear case. 
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Table 3. Comparisons of slope angles, 6,, and deflections, {y,,y,}, calculated from Shifted and Curved 
Displacement Transfer Functions for the P=600 lb nonlinear case; calculated depth factors, c, (second 


column of table 2), and nonlinear lower surface strains, €,, of table 114 of Appendix I were used as 


inputs. 
Shifted Displacement Curved Displacement 
Transfer Functions Transfer Functions 
i ¢, im €,, in/in tan 0, 6,, deg y;, in 6,,rad_ | 0.,deg y,, in 
Calculated Nastran nonlinear Eq. (6a) Eq. (6a) Eq. (6b) Eq. (18a) Eq. (18a) Eq. (18b) 
Eq. (21) lower surface strain 
(table 114) 

0 4.002 0.01155 0.000 0.000 0.000 0.000 0.000 0.000 

1 3.852 0.01167 0.044 2.540 0.330 0.044 2.541 0.330 

2 3.703 0.01177 0.091 5.194 1.342 0.091 5.208 1.342 

3 3.554 0.01185 0.140 7.954 3.069 0.140 8.006 3.069 

4 3.406 0.01191 0.191 10.810 5.545 0.191 10.940 5.545 

5 3.258 0.01195 0.245 13.748 8.809 0.245 14.018 8.809 
6 3.109 0.01194 0.301 16.749 12.898 0.301 17.243 12.898 
7 2.961 0.01190 0.360 19.792 17.851 0.360 20.619 17.851 
8 2.814 0.01182 0.421 22.855 23.708 0.421 24.149 23.708 
9 2.666 0.01168 0.486 25.911 30.509 0.486 27.835 30.509 
10 2.519 0.01148 0.553 28.935 38.296 0.553 31.675 38.296 
11 2.371 0.01120 0.622 31.899 47.107 0.622 35.662 47.107 
12 2.224 0.01084 0.694 34.776 56.980 0.694 39.785 56.980 
13 2.077 0.01038 0.768 37.539 67.949 0.768 44.027 67.949 
14 1.931 0.00980 0.844 40.163 80.040 0.844 48.355 80.040 
15 1.785 0.00909 0.920 42.620 93.271 0.920 52.724 93.271 
16 1.641 0.00819 0.996 44.881 107.643 0.996 57.057 107.643 
17 1.501 0.00706 1.069 46.900 123.132 1.069 61.227 123.132 
18 1.367 0.00559 1.135 48.611 139.668 1.135 65.015 139.668 
19 1.264 0.00355 1.187 49.883 157.104 1.187 68.001 157.104 
20 1.000 0.00037 1.2124 50.464 | 175.138* | 1.2124 | 69.416 175.138* 


*Tdentical Identical 


Note also from table 3 that the deflections, {y,,y,}, calculated respectively from the Shifted 


deflection equation (6b) and Curved deflection equation (18b) using Nastran nonlinear surface strains, 
turned out to be identical. Note also that the Curved slope angle, (0;)..,,,.7 » Which is the true slope angle 


of the deformed embedded beam, reaching up to (0,)..,,,.¢ ~ 69 deg at the beam tip. On the other hand, 
the Shifted slope angle (0,) shined is slightly smaller because the Shifting process reduces the slope angle, 


causing the Shifted slope angle to reach up to only (9,,) shifted = 50 deg at the Shifted beam-tip (fig 5). 
However, as shown in table 3, the slope (tan 6; ) cried calculated from the Shifted slope equation (6a) has 


exactly the same value as the corresponding slope angle (@;),.,..7 Calculated from the Curved slope 
equation (18a) [thatis, (ta 9;) shied = (A Veuved | 
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Figure 16 shows the span-wise plots of slope angles {(9;) snifea (9; )cuveat and deflections {y,,y,} 
based on the data listed in table 3 for P=600 Ib Nastran nonlinear case. Note from figure 16 that the 
deflection curves of {y,,y,} form a single curve because (y, = y,). The (9:) shifted -curve practically falls 
on top of the (;)¢,.q-Curve in the inboard region up to strain-sensing station 7 = 4, and then gradually 


diverge downward from the (@,),.,,,,q -curve toward the beam tip (=n). 


Deflection Identity (y, = y,) 


The unexpected discovery of the deflection identity, y,=y, (table 3, fig. 16), indicates that the 
Shifted vertical deflection, y,, is actually the straightened version of the Curved deflection, y,. This 
deflection identity, y, = y,, can be explained as follows: The Shifting of the deformed material points 
(lying along the neutral axis) to their respective original undeformed x-positions [that is, u > 0] causes 
(dy/dx)’ — 0 (details in Appendix A, ref. 9), thus the Lagrangian curvature equation (4) is reduced to 
the mathematical form similar to the Curved curvature equation (10) as shown below in equation (28): 


Increased to maintain 


same €(s)/c(s) 
1 : d’y/dx? z d’y/dx? _ d’y = d’y _ €(s) (28) 
RQ) fl-(dy/dx? — Tiatantsoy 1-0 dx dk’ es) 


Increased Shifted — Curved 


Equation (28) shows that, based on the same undeformed x-system, the correct Shifting process will 
cause the value of y to match the value of y to maintain the same value of €(x)/c(x), resulting in the 


deflection identity y, = y, (table 3). 


Correct Shifting 


The original horizontal shifting of point A’ to point A” shown in figure 5 was found not to be an 
accurate process of obtaining the true Shifted deflection, y , and, therefore, the shifting process need to be 


modified. 
Figure 17 is the amended figure 5 and graphically shows that the correct Shifting is not the horizontal 
Shifting (fig. 5), but is to move point A’ to point A” in such a way as to bend the curved deflection, y , 


into the equivalent Shifted vertical deflection, y . Thus, the Shifted Displacement Transfer Functions, just 


like the Curved Displacement Transfer Functions, are also applicable to shape predictions of structures 
under large geometric nonlinear deformations provided the value of the Shifted slope 
(tan ®;) sriteal= (9; cunea | iS treated as the true slope angle because the shifting process will reduce the 


true slope angle, ZB’A’C , to a smaller shifted angle, ZB”A”C’ (see fig. 3 and table 3): 


Prediction Errors 


Table 4 lists the beam-tip deflections {y,,y,} of the tapered cantilever tubular beam under all 


loading cases calculated respectively from the Shifted [eq. (6b)] and Curved [eq. (18b)] Displacement 
Transfer Functions, and from Nastran linear and nonlinear analyses. The percent prediction errors listed in 
table 4 were calculated from the error equations (26) and (27) respectively for vertical and curved 
deflection cases. 
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Table 4. Comparisons of beam-tip deflections, {y,,y,}, of the tapered cantilever tubular beam, 


calculated from Shifted and Curved Displacement Transfer Functions and from Nastran linear and 
nonlinear analyses. 


Straight deflection y, , in Curved deflection y, , in 

Shifted Curved 
pe | woe | DaRMERT | ons» | Nagun | Dispacemem | oon 
(Reference) Functions Eq. (26) (Reference) Functions Eq. (27) 

Eq. (6b) Eq. (18b) 
50 23.406 23.424 0.0769* 23.210 23.220 0.0431 
100 46.811 46.828 0.0363 45.337 45.374 0.0816 
200 93.623 93.669 0.0491 83.903 84.082 0.2133 
300 140.434 140.502 0.0484 114.508 114.820 0.2725 
400 187.245 187.314 0.0369 138.776 139.187 0.2962* 
500 234.057 234.161 0.0444 158.439 158.866 0.2695 
600 280.868 280.970 0.0363 174.763 175.138 0.2146 


**Peak errors 


Note from table 4 that at P=600 Ib, the beam-tip curved deflections, y,, (predicted or Nastran- 
generated) is only about 62 percent of the corresponding beam-tip vertical deflections, y, (predicted or 
Nastran-generated). Note also that for vertical deflection cases, the prediction errors are extremely small, 
in the range of 0.0363 percent to 0.0769 percent, with a peak prediction error of 0.0769 percent that 
occurred at P=50 lb. For the curved deflection cases, the prediction errors are slightly larger than the 
vertical deflection cases, but still in the small range of 0.0431 percent to 0.2962 percent, with the peak 
prediction error of 0.2962 percent (3.85 times the peak prediction error of 0.0769 percent for the vertical 
deflection case) occurring at P=400 Ib. 

Figure 18 shows the prediction errors listed in table 4 plotted as functions of applied beam-tip load P 
for vertical and curved deflection cases. For the vertical deflection cases, the prediction error curve is 
almost horizontal, indicating that the prediction error is practically insensitive to the applied load, P. For 
the curved deflection cases, the prediction error increases convex upwardly with applied load, reaching 
the peak of 0.2962 percent at P=400 Ib, and then slightly tapering down to 0.2146 percent at P=600 Ib. 
Keep in mind that the Nastran curved deflection (V),,y appearing in the prediction error equation (27) 
(fig. 5) is calculated from the straight-to curve-deflection conversion equation (20), which gives circular 
arc deflection for (V)yy . However, the theoretical curved deflection y is not a true circular arc. The 


slight decrease in the prediction errors beyond P=400 Ib (fig. 18) could be attributed to the slight decrease 
in the difference between y (non-circular curve) and (¥V),y,y (circular curve) at increasing bending with 


outboard region bend more than the inboard region because of the tapered beam (table 4). 


DISCUSSIONS 


To use the Shifted and Curved Displacement Transfer Functions to calculate out-of-plane deflections 
for structure deformed shape predictions, the following input parameters are required. 
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Lower and upper surface strains, {€,,€,} (= 1,2,3,...,1) —measured at strain-sensing stations 


evenly distributed along the lower and upper strain-sensing lines on the embedded beam surfaces 
(assuming the location of neutral axis of the embedded beam is unknown). For the calculations of 
overall deformed shape of the structure (For example, aircraft wing) under bending and torsion 
loading, an additional embedded beam is needed to form four-line sensing system. 


Domain lengths, Al/(=//n)—once the distribution of strain-sensing stations is defined, the 


domain length, AZ, (strain-sensing station separation distance) is specified. 


Embedded beam depth, h,, at strain-sensing stations i (known for a given structure). 


Depth factors, c, (@=1,2,3,...,.2)—usually unknown for complex structures, and must be 


calculated from equation (21) for each embedded beam using pairs of lower and upper surface 
strains, {€,,€,} , and the embedded beam depth, h,. 


If the depth factors, c;, are known, only the lower surface strain-sensing line is needed (fig. 3). For 


very large geometric nonlinear deformations, the neutral axis can shift with the load level and, therefore, 
the depth factors must be constantly updated using equation (21) at each loading level. 


CONCLUDING REMARKS 


Using the true curved deflection, y, the embedded beam curvature-strain differential equation was 


piecewise integrated to formulate the Curved Displacement Transfer Functions for geometrical nonlinear 
large deformation structure shape predictions. Nastran linear and nonlinear analyses were performed on a 
tapered tubular cantilever beam to analytically obtain surface strains needed for inputs to the Shifted and 
Curved Displacement Transfer Functions for shape predictions. The Nastran-generated deflections were 
then used as a validation reference to study the shape prediction accuracies of the Shifted and Curved 
Displacement Transfer Functions. Some highlights of the results are listed below. 


1. 


For large deformations, one must use the true curved deflection, the curvilinear distance traced by 
a material point from its undeformed position to its deformed position. The traditional vertical 
deflection is merely a vertical component of the true curved deflection. 


The vertical deflections calculated from the Shifted Displacement Transfer Functions and the 
corresponding curved deflections calculated from the Curved Displacement Transfer Functions 
were found to be exactly the same, implying that the vertical deflections based on the Shifted 
formulation are actually the straightened version of the curved deflections based on the Curved 
formulation. 


By replacing {¥,,Qomeat With {y,,(tan@;)s,ic.at» the Curved Displacement Transfer 
Functions can become the Shifted Displacement Transfer Functions. 

Both the Shifted and Curved Displacement Transfer Functions are very accurate, and are 
applicable to the shape predictions of the cantilever tubular beam under geometrical nonlinear 


large deformations with beam-tip deflection reaching as high as 58 percent of the span, for which 
the beam-tip slope angle reaches 69 deg. 
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The Shifted and Curved Displacement Transfer Functions are purely geometric in nature, and 
therefore, one can compute the correct deflections whether the input surface strains come from 
linear or nonlinear deformations. 


For vertical deflection cases, the prediction error of using the Shifted Displacement Transfer 
Functions is in the negligible range of (0.0363—0.0769) percent, and is practically insensitive to 
the change of applied load. 


For the curved deflection cases, the prediction errors of using the Curved Displacement Transfer 
Functions are in the small range of (0.0431—0.2962) percent, slightly larger than the vertical 
deflection cases; because Nastran curved deflection is a circular arc, but the theoretical curved 
deflection is not a circular arc. 


For nonlinear deformations, the magnitudes of the lower and upper surface strains at the same 
strain-sensing station can be slightly different because of the curved-beam effect, which induces 
neutral axis shifting and small axial strain components. 


For nonlinear deformations, either the Neutral-Axis-Shifting Method or the Axial-Strain- 
Elimination Method can be used to obtain correct data (depth factors, bending strains) for input to 
the Displacement Transfer Functions for shape calculations. 


FIGURES 


ED03-0180-02 


Figure 1. A super-long flying wing Helios prototype (wing span 247 ft) under very high wing dihedral 
deformation just before breaking up. 


ED03-0180-03 


Figure 2. Helios Prototype broke-up in mid-air on June 26, 2003 at 2800 ft altitude under very large wing 
dihedral deformation (wing tip deflection reaching 40 ft), and fell into the Pacific Ocean. 
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Figure 3. Embedded beam (depth-wise cross section) of the tapered cantilever tubular beam with strain- 
sensing stations evenly distributed along lower and upper surface strain-sensing lines. 


Incremental 


Beam elastic curve 
beam segment ae 


(deformed neutral axis) 


A'B'= AB[1 + e(s)] - 
Se eas Sea eaae Bie et ag ee a oak 


“~— Undeformed 


OA’ _ 4) AB. 4 45/9) ~ 1 = 89) 


OA R(s) AB R(s)  c(s) 170000 


Figure 4. Small segment of a deformed nonuniform embedded beam for geometrically relating local 
radius of curvature, R(s), to associated surface bending strain, €(s), and depth factor, c(s) . 
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Figure 5. Actual elastic curve of a deformed embedded beam, showing true curved deflections, ) ; the 
curvilinear distances traced by the material points from their undeformed positions to respective deformed 
positions. Horizontal Shifting converts curved deflection, ). into vertical deflection, ¥. with reduced 
slope angle (ref. 15). 
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Figure 6a. Shell-element model for linear analysis. 


Figure 6. Undeformed and deformed shapes of Nastran models of the tapered cantilever tubular beam 
subjected to beam tip load, P, for linear and nonlinear analyses. 
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Figure 6b. Beam-element model for nonlinear analysis. 
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Figure 7. Graphically converting Nastran nonlinear deflection outputs, {u;,y,}, into a circular arc 
length, y,, for comparison with the theoretical curved deflection, y,. 
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Figure 8. Using lower and upper surface strains, {€,,€,}, to obtain lower and upper depth factors, 
{c;,¢,}, based on neutral axis shifting in nonlinear bending. 
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Figure 9. Using lower and upper surface strains, {€,,€,}, to obtain true bending strains, (€,—€,)/2, 
eliminating axial strain components induced in nonlinear bending; known depth factor, c, = h,/2. 
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Linear case P =600 Ib 
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Figure 10. Span-wise distributions of lower surface strains, €,(=—€,), generated from Nastran linear 
analysis of the tapered cantilever tubular beam subjected to different beam-tip load, P. 
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Figure 11. Comparisons of vertical deflection curves for the tapered cantilever tubular beam under a 


different beam-tip load, P, calculated from the Shifted Displacement Transfer Functions [eq. (6b)] and 
calculated from Nastran linear analysis. 
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Figure 12. Span-wise distributions of lower and upper surface strains, €;(2—€,), generated from Nastran 


nonlinear analysis of the tapered cantilever tubular beam subjected to a different beam-tip load, P. 
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Figure 13. Comparisons of curved deflection curves for the tapered cantilever beam at a different beam- 
tip load, P, calculated from the Curved Displacement Transfer Functions [eq. (18b)] and calculated from 


Nastran nonlinear analysis. 
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Figure 14. Plots of known and calculated depth factors, c,, for the tapered cantilever tubular beam 
showing slight neutral axis shifting due geometric nonlinear deformations, especially in the outboard 
region; P=600 Ib nonlinear strain case. 
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Figure 15. Linear (vertical) and nonlinear (curved) beam-tip deflections, {y,,y,}, of the tapered 


cantilever tubular beam plotted as functions of the applied beam-tip load, P. The theoretical deflection 
curves [eqs. (6b) and (18b)] and the corresponding Nastran-generated linear and nonlinear deflection 


curves are graphically indistinguishable. 
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Figure 16. Span-wise plots of slope angle, @,, and deflections, {y,,y,}., of the tapered cantilever 


tubular beam calculated from the Shifted and Curved Displacement Transfer Functions [eqs. (6a), (6b), 
(18a) and, (18b)] using the same Nastran-nonlinear strain data for P=600 Ib listed in table 3. 
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Figure 17. Correct Shifting of deformed material points to their respective undeformed x-locations to 
bend curved deflection, y, into an equivalent straight (vertical) deflection, y (= y); actual slope angle, 
O(s) , turns into slope, tan@(x)[= @(s)], for the Shifted case. 
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Figure 18. Plots of linear (vertical) and nonlinear (curved) beam-tip deflection prediction errors as 
functions of beam-tip load, P. 
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APPENDIX A 
DERIVATIONS OF SLOPE ANGLE AND CURVED DEFLECTION 
EQUATIONS IN RECURSIVE FORMS FOR NONUNIFORM 
EMBEDDED BEAMS 


Appendix A shows the mathematical details of stepwise integrations of the slope equation (15) and 
the deflection equation (17) for the nonuniform embedded beams to obtain the final mathematical forms 
given respectively by nonuniform parts of equations (18a) and (18b). 


Piecewise Linear Strain Representations 
For the piecewise integrations of the slope equation (15) and deflection equation (17), both the depth 
factors and surface strains, {c(x),€(x)}, in the domain x,, <x <x, between the two adjacent strain- 
sensing stations, {x,,,x;}, can be expressed with linear functions given respectively by equations (A1) 
and (A2) [duplications of equations (11) and (12) respectively]: 


e(x)=¢,,-(G4-¢) 4 5 (4%, Sx8x;) (Al) 
Al 
e(x)=€,,-(€.,-€)> a - (x, Sx<x,) (A2) 
Slope Angle Equation 


The slope angle, 0(x), of the nonuniform embedded beam in the domain x,_, <x <x, between the 
two adjacent strain-sensing stations {x,,,x;} is given by equation (A3) [duplication of equation (15)]: 


O(x)= io ade 8, 6, * @4S25%) i 
Slope 


Slope increment at xj 


Substitute equations (Al) and (A2) into equation (A3), and carrying out the integration as follows in 
equation (A4) (ref. 21): 


Xj 
x (E55 E;) Al 
A(x) = f dx+0_, = 
Xj-1 Cc (c, c i= = Xj-1 
i-l i-1 i Al 
Ca(En—€) — EalGa-¢) 
_ a7 E; 476) i ) 


"(Aye ? 
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E, C,- EC, C.-C 
re oa Xj_ prays 8Ei og] Cs Der— a, d+ 6.4] -top. crf +6, (A4) 
i-l i-l i 


At the strain-sensing station, x,, one can write x, —x,, = Al, and equation (A4) yields the slope angle 


6,[= 0(x;)] at the strain-sensing station, x,, as equation (AS): 


E, Cc 
= Fei ayy ap Sei Fifer lo —log,c_,)+94, 
(c.. ~c,) (c,. ~c¢, ie ( Se i Ze 4) i-1 


er 2S. ch) (AS) 


i 


After grouping terms, equation (A5) becomes the final form of the slope-angle equation (A6) for the 
nonuniform embedded beams: 


0, = Al oe oe +6, , 


(c_ 1 —¢,) (¢,. yey Ci 


A 
: G=12,3.4.7) oe 


Equation (A6) is the nonuniform part of equation (18a) in the text. 


Curved Deflection Equations 


The curved deflection, y(x), of the nonuniform embedded beam in the domain x,, <x <x, 
between the two adjacent strain-sensing stations, {x,,,x,;}, is given by equation (A7) [see equation 


17)]: 


1 C(x) (A7) 


y(x)= i dlls EOD a 0. fe iS J O(x)dxt+y¥_, 3 (x, Sx<x,) 
FS 
O(x) 


Substitute equations (A1) and (A2) into equation (A7), and carrying out the integration as follows in 
equation (A8) (ref. 4): 


y(x)= J. J. A +6,, ax+y,= 


“J ET Enie; — EC-1 (¢_1—¢;) ) 
= ee (samme + Al ++ = lo De +¢. = lo pie 4 dx 
fea ; = a ) (Cc, iC; y z.( Al ( 1) i-1 g. i-1 


+ V4 a= Xi_1 )0,4 
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EC) — EC 


_ (ee. 1 ey as ~~ G) 
“Fe.< 5 e Xa +Al ~,, al Al 


Al (A8) 


xf [toe (22 Peon. dhe.s) lose Ca a nid 4, 


At the strain-sensing station, xX, , one can write X,—X,_, = Al, and equation (A8) yields the curved 


(G24) x 


deflection, y,[= y(x,)], at the strain-sensing station, x, , as equation (A9): 


- é. 9 €.,6;—E,C, 
5 = oe =L (Al)? - (Al? Et 5 
AE Ci 1 eo) (¢..- Gy 


=(c.5=¢)) =(C.4=C.) (c,,-—¢;) ss 
x {[-S=Paree., flos. (“Se Bar, }-tos. cn} arta aif. +(Al)6,_, 


S.. e _ (AY) 
= (Al)’ ee eel ¢, (log, G —log.c¢_,)+(c, al} es. +(AL)O,_, 
(F=1,2,35,1) 


Equation (A9) can be written in the final form of the deflection equation for the nonuniform embedded 
beams as equation (A10): 


=; Ej Cj — EiC,_ =| 


Ci _ 
_ ee ie GL ee ADO. 
=(Al)’ Be ey dene: clog, —— + (¢. col}. +(AD6,, 


i-l 


(A10) 


@=1,2,3,...,) 


Equation (A10) is the nonuniform part of equation (18b) in the text. 
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APPENDIX B 
DERIVATIONS OF CURVED DEFLECTION EQUATIONS IN 
SUMMATION FORMS FOR NONUNIFORM EMBEDDED BEAMS 


Appendix B presents mathematical steps to obtain the final dual summation forms given by the 
nonuniform part of equation (18c). The slope-angle equation (A6) and the deflection equation (A10) in 
recursive forms for nonuniform embedded beams are duplicated below as equations (B1) and (B2). 


Slope angle equation: 


E. & =€.¢C. G 
6.=Al aan let Hog, +6, Y= 1,2 Baht) (B1) 
(¢, 1 ¢) (c,, =e) Ci 
Curved deflection equation: 
5 = (ANP ye Sets Et) 6 tog + (c,,-6,) |b + Fe + (ADO, (B2) 
2(¢4-¢;) (G4 —c¢,) i-1 
(i =1,2,3,...,n) 


Equations (B1) and (B2) can be combined into a single deflection equation in dual summation form 
as follows. Writing out equation (B2) for different indices, i, and making use of the indicial relationships 
expressed in equations (B1) and (B2), one obtains equations (B3) through (B6): 

For i= 1: 

é Ege, — E,C bs 
= (Aly a a ea log, — “(a -0) + y, +(AD6, (B3) 
Wea-¢) (G-4) Cy 
For i = 2: 


= (Ap)? 4 £2 FIT) ogg 2 4 (6, -6,) $45, + (ADS, 
Ac, = ra 


C) (G=6) 1 


rT & Eg — ELC 
= (Aly [f= Gy) 84 108. Ge o} 


yi 


(My) Oe Fok i “fc log, +(e, — 7) +5, + (ADO, 
ve -c) (q@-«) Ci 
(Al 8, 
+(ap7| B28. Es og “} scan, 
Co-C, (Cp =a) i. 


40 


E €,C5: E56 CG &—€E En; EC C. 
=(Aly = W232 He log, =+(c ca} ! fo! |. log, ++(c =<) 
{2 CG) (c,-¢,) ° Gj ome). (qy-¢,) ei, os 


Eo ~E, , E0Ci= 


£1 Jog, a Jy + 2(AI)O, 
C=c,- | (=e iy Ce 
(B4) 


E; _ €4C3 Ex 7 
; = A 
= (Al) | 8s fc, =2,) c og, . (eG; cpp 1)6, 


= (Aly a ee £30 [. log, ey (c,- <9] 
Ae =0,). (=a) Cy 


Jo 


(Aly? E, — €, E,Cy — EC, He log, Pm ne co} & — €, EC, — E,Co aE log, 1 F + (Cy -«)| 
0 


4e-¢) (4-4) 1 W(cy-6,) (GQ) 


+ (Aly f= A 4 Foi 7 Fi% tog © Jy + 2(A)0, 
Coe, (Cp —¢)) Co 
(Al), 
+(apy| 28 4 A127 F261) pe pro Ss » Soa E,Cy log, “can, 
C=, e—c,) CG MTG (coe) &o 


€, —€, EC, — E,C, 
- log, 2+(c—-c 
Ac, —¢) (c, —C, y c g, C, ( 1 »| 


9&3; Ey €3 — E5Cy _ 
= (Al) {=e nee fe ; log, ate co} 


EE. E,C; —&,G, 
+—*% _1__ 1 19) ¢ log, ++ (c, -c¢ 
WG =€,) (c,-¢,)° 1 g,—t Cy ( 0 »} 


E,C, - 


$2 081 F160 Jog “fs Jy + (ANB, 
C=6 (eG c ce, (cy -¢,) an 
(B5) 
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a €,-—€, €,C; = EC: Cc = 
M4 =(Al)’ ar 4 = C4 log, ++ (c;—¢,) + y; + (ADO, 
2(cy—C,) (C,—€4) C3 


E,—€E E,C, — E,€ Cc 
A | one (c;-¢,)° ame Es ce »|| 


¥3 


€,=€ 6,65 ='ExC Cc € = & E.C5-= EC Cc 
Aly « —2—"3 -_— -2°'3 "'3"2) © log, 3+ (c, —c,) |+—1 _- 1? "11 6 log, 2 4+(¢ -¢ 
( ) {32 =) (c, -¢,) 3 Se C, ( 2 >| 2(c, ~c,) (c, -¢,) 2 Se Cc ( 1 >| 


1 
E,—E EC, — E,C, c 
+ +—0 1-01 | clog, ++ (cy) -¢,) 
He =c,) (ce) Gy 


& = € E,C, — E4€ 
+(an| 1 Zi 1~2 oe lo 
C) 7 Cy (c,-C,) 


£1 Jog, = if Jy + 3(AN)O, 


1 Co —C; (=e) 0 


=£ E,C, — EC. € é,—€& E,C, — E,€ 
car| 34 203 32 | 3 1 2 1M 2‘1 


2 e + + 2 
C= (oc) &. €=6. (6-6) 


Cc. 
log, — 
C, 


F100 Jog “cane, 
Cy | (cy =G) Co 


w(ar| Sn Be, log.“ 


é, —€ E,C, — E,C c 
5 SL r(eney fro — con, 
1e=c;) (eo. €,) Cc 


’ eG). 4G= C3) 


€ = €& E05. Es. Cc €5 —€ EC; — EG; Cc 
+ —* 4 = 1 8g, log, ~+(4 =e) |+—* = + | g log, — + =e) 
Wo=6) c=c,) es HC. =G,) c 


(c =¢) 


+(Aly? f6 3 4 Bala FH Jog 83 49 F1 Fo 4 9 Fie Flt og 43 Fo Fi 4 3 FoF Fi6o tog “| 
C,—C,  (C) —€3) “Cy Ci) — Cy (c, —¢, GQ %-e (cy —¢;) Co 
+y, + 4(Al)6, 
(B6) 


Based on the indicial progression patterns in equations (B3) through (B6), one can write the curved 
deflection, y,, in a generalized form with two summations (with different summation limits) as equation 


(B7): 
i See Cc, 
=(Al a eles. 4-8 
=( vores ~c;) (c,,-¢,) c Big - 


E 
4(ATy Sei a A A ere oes 
Cj 6; 


j=l (C4= ey 


(B7) 


j-l 


Les Yo + (ADO, 


(i=1,2,3,....n) 
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Equation (B7) is the nonlinear part of equation (18c) in the text. A set of three equations (B1), (B2), and 
(B7) are called Curved Displacement Transfer Functions for the nonuniform embedded beams. 

It must be mentioned that equations (B1), (B2), and (B7)} or eqs. (18a)—(18c) in the text cannot be 
applied directly to the uniform embedded beam because, in the limit of (c,_, =c;), the logarithmic terms 
and the denominators with (c,,—c;) factors will go to zero [that is, log,(c;/c;,_,)=0., (c_,-¢,) =901], 
causing mathematical indeterminacy. Therefore, for the uniform embedded beams, slope and deflection 
equations can be derived separately using a constant depth factor, c(x)=Cc, in the piecewise integrations 


of equation (13) as shown in Appendix C. 
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APPENDIX C 
DERIVATIONS OF SLOPE ANGLE AND CURVED DEFLECTION 
EQUATIONS IN RECURSIVE FORMS FOR UNIFORM 
EMBEDDED BEAMS 


Appendix C presents the details of integrations of the slope angle equation (15) and the curved 
deflection equation (17) for the uniform beams based on piecewise linear strain representations to obtain 
the uniform parts of equations (18a) and (18b) for the uniform embedded beams. 


Slope Angle Equation 


For the uniform embedded beam with constant depth factor, c(x)=c, the slope angle equation (15) 
[or equation (A3) of Appendix A] for the domain x,, <x <x, between the two adjacent strain-sensing 
stations, {X,_,,x,;} , becomes equation (C1): 


= [* 2 3 a Sxex 
O(x) I. : dx + 6; ; a S22x) (C1) 
ue“, 


Based on the linear representation of surface strains, E(x), in the domain x,, <x <x, shown in equation 
(C2): 


&xX) =E_, —(E., — &) —- 5 (44 Sx S%;) (C2) 


equation (C1) can be integrated to yield equation (C3) (ref. 21): 


1 fx =x. 
a(x)=—] le. ~(e.4* c) 4 Jars 6, , 
Cox Al (C3) 
1 (x-x_,) 
=—| €,,(x-x,,)-(€,, -€,)—_—* _ |+0. 
‘| 1 1 1 aAI 1 
At the strain-sensing station, x,, one can write x,—x,, = Al, and equation (C3) yields the slope angle, 
6,[= 0(x;)], at the strain-sensing station, x, , as equation (C4): 
1 (Al) ' 
6.=-| €.,(A)=(8,, =) eG. S:> ADS 123 if) (C4) 
Cc 2Al 


After grouping terms, equation (C4) takes on the final form of the slope angle equation for the uniform 
embedded beam as equation (C5): 


Al 
6= 56 eH +€;)+6., ;  (@=1,2,3....,2) (CS) 
Cc 


Equation (C5) is the uniform part of equation (18a) in the text. 
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Curved Deflection Equations 


For the uniform embedded beam with constant depth factor, c(x)=c, the curved deflection equation 
(17) for the uniform embedded beam in the domain x,_, Sx <x, between the two adjacent strain-sensing 
stations, {x,_,,x,;} , becomes equation (C6): 


y(x) = J. i —dx 7 0. Jas r Vix = J. 0(x) dx 3 Vi 3 (Xj = x = x;) (C6) 


Substitute equation (C2) into equation (C6), and carrying out the integrations as follows in equation (C7) 
(ref. 19): 


a 1 x x _ . = x 
y(x)= =| J le. (€,,—€;) sal Jas dx+y,_,+ i 0,_, dx 
CP XAS XA Xi-1 


Al 
1 px (x Xi ¥ = 
= =i fet Aa (Ea Sa |i +4 (4 = 2,4)0 5 (C7) 
1 (7=23 i (Aa. i 
= 1 i-] =o Ea ree +94 +(%-%,)8.4 


At the strain-sensing station, x;, one can write x,—x,, = Al, and equation (C7) yields the curved 
deflection, y,[= y(x;)] at the strain-sensing station, x,, as equation (C8): 


sills 
a 6Al 


| (Aly € _ 2) Ay 
Cc i-l i-l i 


Jes tne, & ALS 1, 2S isvegl) (C8) 


After grouping terms, equation (C8) takes on the final form of curved deflection equation for the uniform 
embedded beam as equation (C9): 


(2e,,+€,)+59,,+(AD6,, ; (@=1,2,3....,2) (C9) 


Equation (C9) is the uniform part of equation (18b) in the text. 
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APPENDIX D 
DERIVATIONS OF CURVED DEFLECTION EQUATIONS IN 
SUMMATION FORMS FOR UNIFORM EMBEDDED BEAMS 


Appendix D presents mathematical steps to obtain the final summation forms of the curved deflection 
equation given by the uniform part of equations (18c) using the slope angle equation (C5) and the curved 
deflection equation (C9) in recursive forms for uniform embedded beams. Equations (C5) and (C6) are 
duplicated below respectively as equations (D1) and (D2): 


Al 
0. = 5 en +€,)+0,, ; (i = 1,25 Aycoigtt) (D1) 
c 


(2e,,+6,)+3,,+(ADO,, 3 (@=1,2,3,...,2) (D2) 


By combining equations (D1) and (D2), one can obtain a single curved deflection equation for y,, 
which can be written out for different indices, i, in the forms of equations (D3) through (D6): 


Fori=1: 
CAD” - (D3) 
= (2e, +€,)+ (ADO, +5, 
For i= 2 
z Al) = 
= = (2e,+€,)+9,+ (ADO, 
Aly Aly? _ (Al? (D4) 
= (26, +8,) = (24) +2) +(AD6, +5, + (€, + e,)+ (ADO, 
y (Al)@, 
Aly Aly Al) Z 
aK (26, +2,)+ (26, +6,)+ 5 (e, +) +5, +2408, 
For i = 3: 
_ Al 
v= os (2e,+e,)+y, + (ADO, 
Al Aly Aly Aly 2 
a (2e, +€;) ( : (2e,+€,) ( | (22, +4,)+ (€) + €,) +) +2(ADO, (D5) 
(Al) (Aly 
+(e + £2)+ a (€ +€,) + (ADO, 


(Al)@, 
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Al Aly Aly Aly ALY 
= BY 26, +8,)+ 2 (26, +6,)+ 2 (06, +6,)+ De, +6,)+2 2 (6,46) 
+ y, + 3(Al)O, 
For i= 4: 
~ Al)’ = 
oss : (2e, +e,)+y, +(AD6, 
(Al) (Al) (Aly 
(Aly o Ja “Go (Pea tea) + (26, + £9) + se (2e, +€,) 
~ 6 304 2 2 
. Maas, (e,+€,)+2 (€) + €,)+ Jy + 3(ADO, 
Aly’ Aly’ Aly | 
i (e,+6,) + (e,+2,) +S (e, +€,)+ (ADO, 
(Al)@3 
_ (Al? (Al) (Al) (Aly 
== ——(2e,+€,)= = (2e,+€,)+ os (2e, +e, )+ rs (2e, +€,) 


Aly Aly Al) - 
oe (e,+€ 427 (e, 46,13 (€) +€,)+ Jy + 4(ANO, 


Based on the indicial progression patterns in equations (D3) through (D6), one can write the 
deflection, y,, in a generalized form with two summations (with different summation limits) as equation 


(D7): 


~DU-De, FE) +p +(D(ADO, 


=0 for cantilever beams 


__ (Aly < 
j,= = Yae,+e)+4 
jal (D7) 


Contributions from deflection terms Contributions from slope terms 
(i=1,2,3,....n) 


Equation (D7) is the uniform part of equation (18c) in the text. A set of three equations {(D1), (D2), 
and (D7)} are called the Curved Displacement Transfer Functions for the uniform embedded beams. 
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APPENDIX E 
DERIVATIONS OF IMPROVED SLOPE ANGLE AND CURVED 
DEFLECTION EQUATIONS IN RECURSIVE FORMS FOR 
NONUNIFORM EMBEDDED BEAMS 


Appendix E presents the details of integrations of the slope angle equation (15) and the curved 
deflection equation (17) for the nonuniform embedded beams to obtain the final mathematical forms 
given respectively by the nonuniform parts of equations (19a) and (19b). 


Piecewise Nonlinear Strain Representations 


For the piecewise integrations of the slope angle equation (15) and the curved deflection equation 
(17), the depth factors and surface strains, {c(x),€(x)}, in the domain x,, <x <x, between the two 


adjacent strain-sensing stations, {x,,,x,;}, are described respectively with the following linear and 
nonlinear functions [eqs. (E1) and (E2)] [duplications of equations (11) and (13a)} respectively]: 


C(x) =¢,, + (¢,-¢;. > {%.,S22%,) (El) 
3€,_,-4&, +€ 2E. +E 
&(X) =E,7 a PAL i+} (x- Xj i) + “aq ta es <x< X;) (E2) 


Improved Slope Angle Equations 


The slope angle, @(x) of the nonuniform embedded beam in the domain x,_, {x <x, between the 


iA 


two adjacent strain-sensing stations, {x,,,x;}, is given by equation (E3) [duplication of equation (15)]: 


(x) = fc OOD joa: 6, 3 Oia S¥S%) 
p(t) aie (E3) 
Slope increment at xj) 


Substitution of equations (E1) and (E2) into equation (E3) yields equations (E4) through (E8): 


3€,,-4€, + €, QE.+€ 


. ac a Ht (x _ a) Ee = 
A(x) = dx+6._ (E4) 
Xj-1 C.— 1 i-l 
Cyt a (A=aey) 
Let 
oe 3€,, —48, + €,,, (E5) 
2Al 
€.,-2€. +€. 
B= i-l 1 i+] E 
2(Al) Ee) 
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C= E7 
Al (E7) 
€=(@— x5) (E8) 
then equation (E4) takes on the following simplified form shown in equation (E9): 
Ee, + AE+ BE” 
O(x)= [7H gE +6. E9) 
9 + Ce ° 
After carrying out integration of equation (E9), one obtains equation (E10) (ref. 19): 
Ei G4 Sis 
“e loBelG.4 + Cé) + Aa Ava 08. (6.1 + Cé) 
Se ) 
a(é) - Integration of Ist an ‘ (E-4) Integration of 2nd term in eq. (E-4) 4 0, 
Take +CEY -2c,_(c_,+C&) +c, log.(c., + cé)] 
Integration of 3rd term in eq (E-4) 0 
2 
é.. om oa A 1B 2 B 
= o =i cz 8 ce Joe, “r Ca)+ ce — Yaa - Ca) = 2 Gali - CG) 
2 
é, ; 1B B 
- [2-484 1 ig C47 OH + 2c +6., 
(E10) 
= aCe. — ACc,_, + Bc7, |[log,(c,, +C§)—log,c_,|+ a6 
1B B 
+ sal +CEY -c?, | —2 ou [(¢_, + CE)—-c,,]+8,, 
= alc €,,-ACc,_, + Bc?, |[log.(c,, +C&)—log,c,,]+ ra 
1B B 
+ 7 oF eee$ +C&*)-2 ccs +0,, 


- aCe. — ACc,_, + Bc?, |[log,(c,, +C§)— log, c,_,|+ a 


i salen +C& |+6,, 
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In view of the definitions in equations (E5), (E6), (E7), and (E8), equation (E10) takes on the following 
form shown in equation (E11): 


6(x) _ (Al)? aes (¢; = wis y ee (3é,_, -_ 4; oF Eva) (c, = C,1) c_,+ (E., -_ 2€; FE i) 2 x 
(C=C) (Al) 2Al Al 2(Al)” 
¢.—¢... (3e.,-4e,+€.,,) Al 
x} log,| c_, ++" (x- x._,) |- log, c_, |- —=— 1 —— (x- x 
2 i-1 Al ( 3) g. i / 2Al CG, —C,, ( i ) 
1(€., —2€.+€. Aly —¢, 
(2.5 E; : Ein) ) ( ) 7 26. (es _ Hts, ) +4 eames _ is y a 0. 
2 2(Al) (C) Al 
Al > 1 1 5 
=| GalG Guay) + Gb 4g 4E, + 6G Gaeta — 26, FE Ca |X 
(¢-¢,4) 2 2 (E11) 
q 
C,-G_ 
x tos (c. 7 a 7X4 } —log, te 
3e€.,-4e +e. _ —-2€E4+€ -¢ 
= Ew E, +E, (x = X,4)+ 22-26 al 96, (x = Ka pe ae ee a? [4 6. 
2(¢; ~ C4) A(c, _ C5) Al 


I, 
ee ee 7 


i-l — 


In equation (E11), the two terms {J/,,/,} can be simplified through grouping terms as follows in 
equation (E12) and (E13): 
2,1 1 2 
I, = E_(C;, ~ C4) + 79 8 p 4é, tT Ein) Me; Ch y Mei + 7 en a 2€; + Ei, )ej-4 
2,3 | , 2 1 2 
=|(G ea) + 5 — C4 )C.4 + a Cia JE — [ 2(¢, —€,4)€,4 +e; le * 5 — CC + Cj; i 


3 3 1 1 
= G mode Ga + hae - me + seis =(26;-6.4)C. 46:4 9 RFit 
(E12) 


1 
= xl @e, _ CJC En ~ 2(2c; ~ Ci )C 18; + Cte 
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. Wey =4e 7 =.) ae ee ae (C= 0.5) 
te ard Tem 
1 
= ene 1 4é, — Ei MG; 7C_y y= (E, = 2E, +E.) Jez es ~ X;4) 
i So 
4 €,,—2E, +E, A- , 
A(c,—c,_,) Al 
1 
= —|[-3(«, “Gg Ce len = [4(c, = €,4)+2¢,; le, aa [-(c, G4) = Cg le ies = 4) 
eG 4) 
€,,-2€, +€,,, (x-x_,) 
fo—C.) Al 
—2€,4+8,, (x-x,) 
= ——_|(2c._, - 3c, )€._, +2(2c,-c_, é,- x—X, q Sel 20 Fi E13 
2(c,-¢.4) [( il i) i-l ( i i i) iG CE |( 4) A(c,—C,.,) Al ( ) 


Substitutions of equations (E12) and (E13) into equation (E11) yields equation (E14): 


Al = 
B= Fo gpl ei GMb 2G) 4 G6 pat aG tog, fe. jin 
1 
7 2¢,, — 3¢,)E;4 + 2(2¢, - ¢,,); - 
Ue=6.,) [( C4306, (+ 2266 ee CE |X x, 


er Cee (Aa) 
A(c,—c,,) Al 


a 6... 
(42S 2,) (E14) 


At the strain-sensing station x,, one can write x,—x,,= Al, and equation (E14) gives the slope, 
6,[= 0(x;)], at the strain-sensing station, x, , as equation (E15): 


Al 
6, =———, | (2c, - 66,266.) FCC lo 
i Ac, = ae ea [( Ci ah i i-l ) i Ci iE E.| gt ome 


Al 


+—— | (2¢,_, — 3¢,€,, +2(2¢, -— 6, )€; - 6&4, ]+ if ta (E15) 
(6-63) 


C; Ein 


I; 


(£=1,2,3,.0.;7) 


The term /, in equation (E15) can be simplified as follows in equation (E16): 
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Ie te, - 3¢,)8,, + 2(2e, -¢,,)8, — C6, ] + Pe Bt 
2(c,-€ 4)" A(c,-c,,) 


Al 1 
= rama ke — 364 5G — C4 fe +[2(2¢,-¢.,)-(¢ -¢.)]&; 


1 (E16) 
« ~ 5 ~ cd ]ef 


Al 
= “ean - 3c, )é€_, -2(3¢,-¢,_)é, +(e, +¢, Eta | 
I i-l 


In view of equation (E16), equation (E15) takes on the final form as equation (E17): 


Al 
0, = Renee Cis iNC,é, 1 =2¢,. 1€; )+¢¢ i-1€ é,,, log, C., 
Al 


7 i Cees a y 


(E17) 


[(Sc, — 3¢,_,)€; — 2(3c, — ¢_é, +(C, + 6)Ein | +9, 


i+] 


(i=1,2,3,....n) 


Equation (E17) is the nonlinear part of equation (19a) in the text. 


Improved Curved Deflection Equations 


The curved deflection, y(x), of the nonuniform embedded beam in the domain x,, Sx <x, between 
the two adjacent strain-sensing stations, {x,_,,x,}, can be expressed as equation (E18) [see eq. (17)]: 


x €(x) " x - 
y(x)= J Ie ee ak fori = J, O@)ar+5., SAK eS) 


(x) 


Substitution of equation (E14) into equation (E18) yields equation (E19): 


as Al 
y(x)= a= c4) [(2c, C4 MCE. — 26, ,€;) + C,C,_ fg \% 


x |" oe.(c.¥8 ae (x 5.) log. c., fas 


¥ [2e,, — 3¢,)e,-4 + 22¢, = 6.8 - Ca] [> Ce- x. )dx 


(E19) 


+ 
2(C, — C.-4) 
—2€, Fe ft (x- 

re 7 o_ ) 


x;- ide 
rv dx fo G4 det Yi 


Carrying out the integration of equation (E19), one obtains equation (E20) (ref. 21): 
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30) = — —{le«- C4 MCE — 2C,4€)) + CCE: |X 


Al CG C=G 
x c_, t+ (x-x,) log} c_, ++ (x- x, 
| -C., ( i-1 Al ( ») 2.( i-l Al ( ») 


a “ee -x,)-¢., loge.) (x—x,,)log,c,, } 


(E20) 
1 
+ 2 
A(c, —¢,.) 
6.4 ~ 26 FE (A= x) 
12(c,-¢,,) Al 


[ens = 36 )E 9 t226 =¢., Jb —Ce.. | (x-x,)° 


+(x-%,)O., + Yi 
6 re ey 


At the strain-sensing station, x, , one can write x, —x,, = Al, and equation (E20) gives the curved 


deflection, y,[= y(x,)], at the strain-sensing station, x, , as follows in equation (E21): 


é (Al)’ 
ar ere eee Ce, 1 2G 1€; J+GC i- iE E,,, |X 
x [<, log, c, —(¢, —¢_,)—¢,, log, ¢,, — (¢, - he Cy |} 
(Aly 

to | (264 — 36; JE + 220; — 6 JE, — CE + +6. 1 2E, +E, )(C,—C.,) 
A(c,-¢,_,)° 

t+y,_,+(AD@,, 

(Aly 
Ae. 4 : 


2 
+ =e (26, ~364<i= “ste, +( 4 —2¢.,— ee +s Je +(-« i. = Je 
= as 5° 3 3 3 


+ ¥ P(ADG., 


[(2c, 7C.y MCE ~ 26,4; ) HCC) E 541 | [c, (log, Ci - log, Ci ) = (¢; mol )| 


@=1,2,3,....n) (E21) 


After grouping terms, equation (E21) becomes equation (E22): 


7 (Aly C, 
i eae" ~~ Giles, 4 — 26,56) eee | C; log, C i Cia) 
A 2 
= 2 AAW Tipe. — 50, )E4 —2(5¢, — 20, )€; + (2c, + Ens | +5. + (ADO, (E22) 


12(c,-¢_,)° 
(= 1,253,428) 


Equation (E22) is the nonuniform part of equation (19b) in the text. 
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APPENDIX F 
DERIVATIONS OF IMPROVED CURVED DEFLECTION EQUATIONS 
IN SUMMATION FORMS FOR NONUNIFORM EMBEDDED BEAMS 
Appendix F presents the mathematical derivations of the nonuniform parts of the summation curved 


deflection equation (19c). The Improved slope angle equation (E17) and Improved curved deflection 
equation (E22) are duplicated below as equations (F1) and (F2) respectively: 


Al Cc, 
0, = ————__| (2c. -—c._, (c,€,, —2c¢,,€,)+cc,_,€.,, |lo : 
i 2(c, _ om , [( i i-l Xe i~i-l i-l ) i~i-l il &. Cu 
Al 
= Heo, =, 3,4 Ens ~ 2(3¢; 7 Cy E; + (¢; +. Ein | = 6. —_ 
(i =1,2,3,....2) 
Z (Al) C, 
/= ———| (2c. -c._, )(c,€._, —2¢._,€,)+¢¢,_,€, log,——(c,-c, 
Be Ac, - Ga) [( i i-l ( i~i-1 i-l A) C; i-l a C; ge C4 ( i C4) 
(Al)’ o (F2) 
[(8c, -—Sc,J€., —2(5¢, — 2¢,_,)é, + (2c, + 6, Ej, |+ 9. + (ADO, 


~ 12(¢,-¢.,) 
(i =1,2,3,....7) 


Equations (F1) and (F2) can be combined into a single deflection equation in dual summation form as 
follows. Writing out equation (F2) for different indices, i, and making use of the indicial relationships 
expressed in equations (F1) and (F2), one obtains equations (F3) and (F4): 


For i =1: 
a (Aly* Cc 
y, = ena" = 0, JO Sp 226 — 6, \Gye, cel log, G =(c =o, ) es 
(Al) 
~—— | (8e, — 5cy Ey — 2(Sc, — 2¢y JE, + (2c, + Cy Es | + Hy + (ADO, 
Ie, = cy) 
For i = 2: 
= Aly” Cc 
= a z[ (2c, -— ¢ )er€, — 2(2c, — ¢, cs + €5C,€3 |] c log, +-(c, -¢,) 
2\¢;=¢,) rae 
Aly 7 
Se ae =[(8c, —5¢,)€, — 2(Sc, — 2c, )é, + (2c, +0, )€ | + 9, + (ADO, 
I2e,=c,) 
Al) Cc 
= a1. —C,)e,€, — 2(2c, —c,)c,€, + Opciey| c, log, +-(c, —¢,) 
Ac, = C,) Cc 
(Al) 
= ie, eye — 5c, )€, — 2(5c, —2¢, )€, + (2c, +¢,)€5 | 
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Aly’ Cc 
64, = Cy )C\Eq — 2(2, — Cy CoE, + Gre, | c, log, +-(c,-¢) 
2(C, ~~ Co ) Co 

* 2: 
Al = 
ee =[(8c, — Sc, )Ey — 2(Sce, — 2ey JE, + (2, + Cy Es | + Hy + (ADO, 
126, = 64) 

A 2 
ANY J026, ~ Cy JO Ey — (2c, = Cy )EyE, +. CCE |Log, “1 
2(C,- Cy) Co 

Aly 
pa ee = 3¢y Eq — ABC, = Cy Ey + (Cy + Cy Es] + (ADB, 
A4(c, -—Cy) 


(ADO, (F4) 


After grouping terms, equation (F4) becomes equations (F5) and (F6): 


Aly ; 
_ AN 06, ~¢,)C,€, —2(2c, — ¢, CE, + €,€€, |] cy log, +-(c, -¢,) 
2(c, —¢,) ; 
yo = (any 
“De ay Oe — 5¢,)& — 2(5¢, — 2¢,)€, + (2c, + ¢,)E | 
Aly C 
1% ~ Cy )CEq — 2(2e, = Cy Op, +. CCE> |] C, log. +-(C, - Cp) 
i NAG =, ) a 
(Al)’ ; - 
“THe, -a 78 = 30, 66 = 2096, = 2G, Jb (26, Fee, | +5y (FS) 
2 
Seale = Cy CE — 2(2e, = Cy EE, +. CCpE> ]log, . 
ey cary ° ++ 2(Al)O, 
“Rea _ SC. JE 4 2(3c, = Gj + ic + CED | 
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For [=3: 


. Aly 
po) Op asa -9Ge eee, need, || tor (e=65 
Ais = Co) se 
2 
_ _ AN" Ie, —5c,)€, — (Sc, — 2c, é, + (2c, +c, )é, ]+ 7, + (ADO, 
ie.=c,) 
Aly’ 
Efe, = Cy )€y£4 — 2(2Cy = Cy)C Ey + C4C £4 |]; log, S = (C, - cy) 
: 2(c3 -C,) 65 
Ay 
“Ke naples ~5c,)€, — 2(5e, —2¢,)€, + (2c, + C,)€4] 
Aly Cc 
ta 16 —€,)C€, —2(2c, —¢,)c€, + CC,€3 | c, log, +-(c, -¢) 
AG =o) Ci 
2 
=o NED Tee, —5c,)é, — 2(5c, — 2c, )é, + (2c, + ¢,)é; | 
12(c, -¢) ; 
Al)’ Cc 
OD \(2¢, = Cy )CEq — 2(2e, — Cy CoE, + C1CE> Ie log, +-(c,-Cy | 
AAC — 5) Cc; 
+ ‘ F6 
= A ee ~5€y )Ey — 2(Sc, — 2c )E, + (2c, +5 Es ]+ Fp om 
eG) 
Aly Cc 
(Ad) ; [(2c, = Cy CE — 2(2€, — Cy CoE; + CCoE> Jlog na 
Ps (ete oe) ra 
+(2=1) “aay 
~~ (5¢, — 3¢y Ey — 23 ¢, — Cy JE, H(E, + Cy Ed | 
A(c, — Co) 
+2(Al)6, 
‘ 
2 1G, =O, )O,E, — 2026, — CGE, +0, C,€, log. a 
we =C,) GC; 
(Aly 
He, al — 3c, )€, — 2(3c, —¢,)€, H(c, +¢,)€3|+ 
2 
A Ae =i Og = 2126; =C, Ce FE CiCe, Jlog, ail 
AG oi) GC 
a — 3¢y Ey — 2(3c, = Cy JE, + (C, + Cy Es | + (ADO, 
AC = ¢, ) 7 


(AD, 


After grouping terms, equation (F6) becomes equation (F7): 


56 


A 2 
AN (26, —C,)C,E, — 2(2c, — c, )e,€, + C68, | c; log, a (c,—C,) 
_ }2(e;-¢) Cy 
A 2 
“Telit —5c,)€, (5c, - 2c, )€, + (2c, + 64),] 
3 2 
Al)’ 
—~ lee, = €,)C,€ — 2(2cy — C, )c,£, + €3€\€; Jc, log, ee (c,- <) 
“ (C5= Cj) Cc 
(Al)’ 
“Die, aye — 5c, )€, — 2(5c, -— 2c, )é, + (2c, + ¢,)€; | 
Aly’ c 
_AY _T02¢, = Cy Eq — 2(2€, — Cy OpEy + CoE, | c log, +-(c,-c) 
Z\c =,) Ca = FT 
+ “a + Yo (F7) 
————\, | (8c, — 5c, )é, —2(Sc, — 2c, JE, + (2e, + E> 
Decay [(8q — 509 )€y — 2(Se, —2ey)e, + (2q, +6 )€, | 
Al)’ Cs 
AN" fee, =6))098) = 226; =C))6 85 + CC; |log, —_ 
(3-2) 2(c, -C,) cc 
oe ie — 3c, )é, —2(3c, - ¢,)é, + (Cc, +¢,)€5 | 
A(c, = c,) 2 11: 2, 172 2 17-3) 
(Al)’ Cc 
Henegee” = Cy Ep — 2(2c, = Cy )OpE, + Cp €> flog, a 


+(3=1) + 3(Al)0, 


ee ee ee cer Cee 
A(c, _ Cy y 1 0 0 1 0 1 1 0 2 
Observing the indicial behavior, equation (F7) can be generalized for index i, and the deflection, y,, can 


be expressed in a generalized form with two summations (with different summation limits) as 
equation (F8): 


(Al) o 
; en 2,y ee = Cp MC E)4—2¢,48;)+ tien, | clog, +-(c,-c¢,.) 
j _ Cj Ci Ci 7 j 
a 0 
jal (Al) 


remand ka = Se, ,)€ 4 —2(Se; —2c,,)€; + 2c; +. )E pu | 
J jl 


Contribution from deflection terms 


(Aly C, (F8) 
Pe Xe nei =O) CE — 20Gb) Peek |log. te 
ee Ci ; 
*2G-D (Aly + (1)(Al)O, 
j=1 
j “Fea — 36 4JE 4 -2GC; Cy JE, + (C, + Gao 


Contributions from slope terms 


(i=1,2,3,....n) 
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Equation (F8) is the nonuniform part of equation (19c) in the text. A set of three equations {(F1), (F2), 
and (F8)} is then called the Improved Curved Displacement Transfer Functions for the nonuniform 
embedded beams. 
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APPENDIX G 
DERIVATIONS OF IMPROVED SLOPE ANGLE AND CURVED 
DEFLECTION EQUATIONS IN RECURSIVE FORMS FOR 
UNIFORM BEAMS 


Appendix G presents the details of integrations of the slope angle equation (15) and the curved 
deflection equation (17) for the uniform beams based on piecewise nonlinear strain representations to 
obtain uniform parts of equations (19a) and (19b). 


Slope Equation 


For the uniform embedded beam with constant depth factor, c(x)=c, the slope angle equation (15) 
for the domain x,, XS x, between the two adjacent strain-sensing stations, {x,,,x,}, becomes 


equation (G1): 


A(x) = | ——de+ 0.) 3 (4 S"Sx) an 
Xj-1 (Gi 


The nonlinear representation of strain, €(x), in the domain x,, {xX <x, between the two adjacent 


strain-sensing stations, {x,_,,x,}, described by equation (13a), is duplicated below as equation (G2): 


See ey a GY 
2Al 2(Al) a 


eg S25a,) 


i-L 


E(x) =E,,- 


In view of equation (G2), equation (G1) can be integrated to yield equation (G3) (ref. 19): 


1 pe 3€,,- 46, +€, bi—2e.78 
O(x) = -| €,,-— 4 1 x) ¢ Sr? [de 8, 
Cer 2Al 2(Al) 
1 3€,,-4€, + €, bj =26.76 (G3) 
ea €,,(4—-x,_,)-— Hr +H (xx, 1) +0., 
c 4Al 6(Al) 
(tg SPS) 
At the strain-sensing station, x,, one can write x,—.x,, = Al, and equation (G3) yields the slope angle, 
6.[= O(x,; )], at the strain-sensing station, x, , as equation (G4): 
Al 3€,,-—4€, + €, 6, —28, +, : 
6, = al a 7 bp et - _ +6. & G12 Bult) (G4) 
Cc 


After grouping terms, equation (G4) takes on the following final form as equation (GS): 
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Al 
0, a 008 + 8€; —€,,)+ 8, 6G S123 5,8) (G5) 
Cc 


Equation (G5) is the uniform part of equation (19a) in the text. 


Curved Deflection Equations 


For the uniform embedded beam, [c(x)=c], the curved deflection, y(x), in the small domain 


X,, £X <x; between the two adjacent strain-sensing stations, {x,_,,x;}, can be obtained by integrating 
the slope equation (G3) [see equation (C6)] as equation (G6): 


A x x E(X _ x - 
y(x) = I. i EX) 5. + 0. Jos +y,= I. O(x) dx+ yi : (tes <x< X;) (G6) 


O(x) 


Substituting equation (G3) into equation (G6), and carrying out integrations, one obtains equation 
(G7) (ref. 21): 


= 1 x 3€,, —4€, +, 1 2 E) —2€, rE; 1 3 
x)=- €_,(x-x_,)-— W(x - x, + (x- x, dx 
y( ) ‘| aC 1) A4Al ( i) 6(Al)’ ( i) 
+f 0,_,dx + Yi. 
st | (G7) 
lle, 3e€.,-4€. +2 _,-2€. +€ 
ror Hel og SP Ei Ej Cis (x—x,,)) + Sei * Fa yyy 
c| 2 12Al 24(Al) 
+(x—x,,)0,,+ y4 
(4 S2=2,) 


At the strain-sensing station, x,, one can write x,—x,,= Al, and equation (G7) gives the curved 
deflection, y,[= y(x,)] , at the strain-sensing station, x,, as equation (G8): 


i 


ah ) [Sin as Fl SEE | 5 CANA, 


c 2 12 24 (G8) 
(S12, 3,551) 
After grouping terms, equation (G8) takes on the final form as equation (G9): 
_ (Aly = : : 
y,= rc (TE, FOE, 6,4 FY a PADD. . UH 23.008) (G9) 


Equation (G9) is the uniform part of equation (19b) in the text. 
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APPENDIX H 
DERIVATIONS OF THE IMPROVED CURVED DEFLECTION 
EQUATION IN SUMMATION FORM FOR UNIFORM EMBEDDED 
BEAMS 
Appendix H is to derive the final summation form of the curved deflection equation for the uniform 


embedded beams using nonlinear strain representations. Equations (G5) and (G9) are duplicated below as 
equations (H1) and (H2), respectively. 


Al 
6, =—(5¢,, +82, -€,,)+ 6, 5 @=1,2,3,...,n) a”) 
12c 
_ _ (Al) 
a IAc —— (7, + 6€, - €,)+ 94 +(ADO, 5 @=1,2,3,...,0) af 


Equations (H1) and (H2) can be combined into a single deflection equation in summation forms as 
follows. Writing out equation (H2) for different indices, i, and making use of the indicial relationships 
expressed in equations (H1) and (H2), one obtains equations (H3) and (H4): 


Forj=1: 
j= (Al) 5, + (AN)O, (H3) 
24c 
Fori = 2: 
y, = ay wie +66, —£,)+ 9, 4+(AD@, 
a (7é,4 62, —€, ay Vy + (ALO, (H4) 
ao “De +86, —€,)+(AD6, 
(ADO, 
After grouping terms, equation (H4) becomes equations (H5) and (H6): 
= (Aly (Al) (Al) 
y, = Ade -— (Je, P0f; >) + ihe (7&, + GE, = €,)4(2—1) oy (06,4 8e,~€,) (HS) 


+ y, + 2(Al)6, 
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Fori = 3: 


2 
y3 = a + 6€,—€,)+ y, +(ADO, = oe (7€, + 6€,—€,)+ 
Cc 
Al Al Al a 
ore 2 S Zi “ sy S€5)F 5 t ZAG, 
H6 
a S _ 6.) +(AD8, (H6) 
(AD ®, 


After grouping terms, equation (H6) becomes equation (H7): 


2 2 
y3= ay (7é, + 6€, Ey) = e,- e+ (7€) + 6€, —E,) 
(ALy’ (Aly’ ea 
+ (3-2) (56, + BE — £3) + (3-5 + 86, — €) + Jo + 3(AN)O, 
Cc Cc 


Observing the indicial behavior, equation (H7) can be generalized for index, i, and the curved 
deflection, y,, can be expressed in a generalized form with two summations (with different summation 


limits) as equation (H8): 


5 _ (Aly 
yi “Fe 


2 i-l 
ca ao [SE )4 + 8E; — Ej.) + DADO, 


SEL, FOE, -E,,)+ y+ 
(H8) 


Contributions from deflection terms Contributions from slope terms 


(i=1,2,3,....n) 


which is the uniform part of equation (19c) in the text. A set of three equations {(H1), (H2), and (H8)} are 
called the Improved Curved Displacement Transfer Functions for the uniform embedded beams. 
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APPENDIX I 
SUMMARY DATA OF LINEAR AND NONLINEAR DEFORMATION 
ANALYSES OF THE TAPERED CANTILEVER TUBULAR BEAM 


The complete set of the strain and deflection data generated by Nastran linear and nonlinear analysis 
of the tapered cantilever tubular beam are listed in tables I1—I14 for different values of beam-tip load, P. 
Tables I1-I7 are for linear cases, and tables I-8—I14 are for the nonlinear cases. In the last columns of 
tables I1—I14, theoretically predicted deflection data are listed for comparisons. 


Linear Analysis 


In the Nastran linear analysis, the outputs give only vertical deflections, y,, and zero axial 
displacements (that is, u; =0) (fig 6a). Note from tables I1—I7 that for linear cases, the lower and the 


upper surface strains at the same strain-sensing cross section have the same magnitudes. In the last 
columns of tables I1—I7, the data of Shifted (vertical) deflection, y,, were calculated from the Shifted 


Displacement Transfer Functions [eq. (6b)] using the known depth factors, c,(= h,/2), and lower surface 
strains, €,, listed in tables I1—I7. 


Table I1. Vertical deflections, y,, of the tapered cantilever tubular beam calculated from Nastran linear 
analysis and from Shifted Displacement Transfer Functions [eq. (6b)] using the known depth factors, 


c,(=h, : 2), and Nastran linear lower surface strains, €, ; P=50 Ib at the beam tip. 


€,, in/in €,, in/in y,,in y,,in 
1 c; ,in Nastran-linear upper Nastran-linear lower Nastran-linear Shifted DTF 
surface strain surface strain vertical deflection Eq. (6b) 
0 4.0000 -0.00124 0.00124 0.000 0.000 
1 3.8500 -0.00127 0.00127 0.036 0.036 
2 3.7000 -0.00130 0.00130 0.146 0.145 
3 3.5500 -0.00134 0.00134 0.335 0.335 
4 3.4000 -0.00137 0.00137 0.609 0.608 
5 3.2500 -0.00141 0.00141 0.973 0.973 
6 3.1000 -0.00144 0.00144 1.435 1.435 
7 2.9500 -0.00148 0.00148 2.002 2.002 
8 2.8000 -0.00152 0.00152 2.681 2.682 
9 2.6500 -0.00155 0.00155 3.482 3.484 
10 2.5000 -0.00158 0.00158 4.415 4.418 
11 2.3500 -0.00161 0.00161 5.490 5.494 
12 2.2000 -0.00164 0.00164 6.719 6.724 
13 2.0500 -0.00165 0.00165 8.116 8.122 
14 1.9000 -0.00165 0.00165 9.693 9.701 
15 1.7500 -0.00162 0.00162 11.464 11.475 
16 1.6000 -0.00155 0.00155 13.443 13.457 
17 1.4500 -0.00141 0.00141 15.637 15.655 
18 1.3000 -0.00117 0.00117 18.048 18.069 
19 1.1500 -0.00075 0.00075 20.658 20.679 
20 1.0000 0.00000 0.00000 23 406* 23.424* 


*Extremely close 
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Table I2. Vertical deflections, y,, of the tapered cantilever tubular beam calculated from Nastran linear 
analysis and from Shifted Displacement Transfer Functions [eq. (6b)] using the known depth factors, 


c,(=h, / 2), and Nastran linear lower surface strains, €,; P=100 lb at the beam tip. 


€, , in/in €,, in/in y; ,in y, ,in 
u c; ,in Nastran-linear upper Nastran-linear lower Nastran-linear Shifted DTF 
surface strain surface strain vertical deflection Eq. (6b) 
0 4.0000 -0.00248 0.00248 0.000 0.000 
1 3.8500 -0.00254 0.00254 0.072 0.071 
2 3.7000 -0.00260 0.00260 0.292 0.291 
3 3.5500 -0.00267 0.00267 0.670 0.669 
4 3.4000 -0.00274 0.00274 1.218 1.216 
5 3.2500 -0.00281 0.00281 1.947 1.945 
6 3.1000 -0.00289 0.00289 2.871 2.868 
7 2.9500 -0.00296 0.00296 4.004 4.002 
8 2.8000 -0.00303 0.00303 5.362 5.361 
9 2.6500 -0.00310 0.00310 6.964 6.964 
10 2.5000 -0.00317 0.00317 8.830 8.830 
11 2.3500 -0.00323 0.00323 10.980 10.981 
12 2.2000 -0.00327 0.00327 13.439 13.442 
13 2.0500 -0.00330 0.00330 16.232 16.238 
14 1.9000 -0.00329 0.00329 19.386 19.395 
15 1.7500 -0.00323 0.00323 22.929 22.941 
16 1.6000 -0.00310 0.00310 26.886 26.901 
17 1.4500 -0.00283 0.00283 31.274 31.294 
18 1.3000 -0.00234 0.00234 36.097 36.120 
19 1.1500 -0.00150 0.00150 41.316 41.340 
20 1.0000 0.00000 0.00000 46.811* 46.828* 


*Extremely close 
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Table I3. Vertical deflections, y,, of the tapered cantilever tubular beam calculated from Nastran linear 


analysis and from Shifted Displacement Transfer Functions [eq. (6b)] using the known depth factors, 
c,(=h, 7 2), and Nastran linear lower surface strains, €, ; P=200 lb at the beam tip. 


€, , in/in €,, in/in Y;5 Mi. y, ,in 
u c; ,in Nastran-linear upper Nastran-linear lower Nastran-linear Shifted DTF 
surface strain surface strain vertical deflection Eq. (6b) 
0 4.0000 -0.00495 0.00495 0.000 0.000 
1 3.8500 -0.00508 0.00508 0.143 0.142 
2 3.7000 -0.00521 0.00521 0.584 0.581 
3 3.5500 -0.00534 0.00534 1.340 1.338 
4 3.4000 -0.00548 0.00548 2.436 2.433 
5 3.2500 -0.00563 0.00563 3.894 3.890 
6 3.1000 -0.00577 0.00577 5.741 5.738 
7 2.9500 -0.00592 0.00592 8.007 8.005 
8 2.8000 -0.00606 0.00606 10.725 10.724 
9 2.6500 -0.00620 0.00620 13.929 13.930 
10 2.5000 -0.00634 0.00634 17.659 17.662 
11 2.3500 -0.00646 0.00646 21.960 21.966 
12 2.2000 -0.00655 0.00655 26.877 26.888 
13 2.0500 -0.00660 0.00660 32.463 32.480 
14 1.9000 -0.00658 0.00658 38.772 38.795 
15 1.7500 -0.00647 0.00647 45.858 45.889 
16 1.6000 -0.00619 0.00619 53.771 53.811 
17 1.4500 -0.00565 0.00565 62.549 62.599 
18 1.3000 -0.00469 0.00469 72.193 72.252 
19 1.1500 -0.00300 0.00300 82.631 82.692 
20 1.0000 0.00000 0.00000 93 .623* 93.669* 


* Extremely close 
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Table 14. Vertical deflections, y,, of the tapered cantilever tubular beam calculated from Nastran linear 
analysis and from Shifted Displacement Transfer Functions [eq. (6b)] using the known depth factors, 


c,(=h, 7 2), and Nastran linear lower surface strains, €, ; P=300 Ib at the beam tip. 


€, , in/in €,, in/in Y;5 Mi. y, ,in 
U c; , Nastran-linear upper Nastran-linear lower Nastran-linear Shifted DTF 
surface strain surface strain vertical deflection Eq. (6b) 
0 4.0000 -0.00743 0.00743 0.000 0.000 
1 3.8500 -0.00762 0.00762 0.215 0.213 
2 3.7000 -0.00781 0.00781 0.875 0.872 
3 3.5500 -0.00801 0.00801 2.011 2.007 
4 3.4000 -0.00822 0.00822 3.654 3.649 
5 3.2500 -0.00844 0.00844 5.841 5.835 
6 3.1000 -0.00866 0.00866 8.612 8.607 
7 2.9500 -0.00888 0.00888 12.011 12.007 
8 2.8000 -0.00909 0.00909 16.087 16.085 
9 2.6500 -0.00931 0.00931 20.893 20.893 
10 2.5000 -0.00951 0.00951 26.489 26.493 
11 2.3500 -0.00968 0.00968 32.940 32.949 
12 2.2000 -0.00982 0.00982 40.316 40.332 
13 2.0500 -0.00990 0.00990 48.695 48.719 
14 1.9000 -0.00988 0.00988 58.158 58.191 
15 1.7500 -0.00970 0.00970 68.787 68.832 
16 1.6000 -0.00929 0.00929 80.657 80.716 
17 1.4500 -0.00848 0.00848 93.823 93.897 
18 1.3000 -0.00703 0.00703 108.290 108.378 
19 1.1500 -0.00449 0.00449 123.947 124.039 
20 1.0000 0.00000 0.00000 140.434* 140.502* 


* Extremely close 
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Table I5. Vertical deflections, y,, of the tapered cantilever tubular beam calculated from Nastran linear 


analysis and from Shifted Displacement Transfer Functions [eq. (6b)] using the known depth factors, 
c,(=h, 7 2), and Nastran linear lower surface strains, €, ; P=400 Ib at the beam tip. 


€, , in/in €; , in/in Y;5 Mi. 74 ul 
U c; ,in Nastran-linear upper Nastran-linear lower Nastran-linear Shifted DTF 
surface strain surface strain vertical deflection Eq. (6b) 
0 4.0000 -0.00990 0.00990 0.000 0.000 
1 3.8500 -0.01015 0.01015 0.287 0.284 
2 3.7000 -0.01042 0.01042 1.167 1.162 
3 3.5500 -0.01069 0.01069 2.681 2.674 
4 3.4000 -0.01096 0.01096 4.872 4.864 
5 3.2500 -0.01125 0.01125 7.788 7.780 
6 3.1000 -0.01154 0.01154 11.483 11.475 
7 2.9500 -0.01183 0.01183 16.015 16.008 
8 2.8000 -0.01213 0.01213 21.449 21.444 
9 2.6500 -0.01241 0.01241 27.858 27.855 
10 2.5000 -0.01268 0.01268 35.319 35.320 
11 2.3500 -0.01291 0.01291 43.919 43.927 
12 2.2000 -0.01309 0.01309 53.755 53.771 
13 2.0500 -0.01320 0.01320 64.927 64.953 
14 1.9000 -0.01317 0.01317 77.544 77.582 
15 1.7500 -0.01293 0.01293 91.715 91.769 
16 1.6000 -0.01238 0.01238 107.542 107.613 
17 1.4500 -0.01130 0.01130 125.097 125.185 
18 1.3000 -0.00937 0.00937 144.387 144.489 
19 1.1500 -0.00599 0.00599 165.263 165.367 
20 1.0000 0.00000 0.00000 187.245* 187.314* 


* Extremely close 
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Table 16. Vertical deflections, y,, of the tapered cantilever tubular beam calculated from Nastran linear 
analysis and from Shifted Displacement Transfer Functions [eq. (6b)] using the known depth factors, 


c,(=h, re 2), and Nastran linear lower surface strains, €, ; P=500 Ib at the beam tip. 


€, , in/in €; , in/in Y;5 Mi. 74 ul 
U c; ,in Nastran-linear upper Nastran-linear lower Nastran-linear Shifted DTF 
surface strain surface strain vertical deflection Eq. (6b) 
0 4.0000 -0.01238 0.01238 0.000 0.000 
1 3.8500 -0.01269 0.01269 0.359 0.356 
2 3.7000 -0.01302 0.01302 1.459 1.453 
3 3.5500 -0.01336 0.01336 3.351 3.343 
4 3.4000 -0.01371 0.01371 6.089 6.080 
5 3.2500 -0.01406 0.01406 9.735 9.725 
6 3.1000 -0.01443 0.01443 14.353 14.344 
7 2.9500 -0.01479 0.01479 20.019 20.011 
8 2.8000 -0.01516 0.01516 26.812 26.807 
9 2.6500 -0.01551 0.01551 34.822 34.822 
10 2.5000 -0.01584 0.01584 44.148 44.154 
11 2.3500 -0.01614 0.01614 54.899 54.913 
12 2.2000 -0.01637 0.01637 67.193 67.217 
13 2.0500 -0.01650 0.01650 81.158 81.195 
14 1.9000 -0.01646 0.01646 96.930 96.982 
15 1.7500 -0.01617 0.01617 114.644 114.716 
16 1.6000 -0.01548 0.01548 134.428 134.522 
17 1.4500 -0.01413 0.01413 156.372 156.489 
18 1.3000 -0.01172 0.01172 180.483 180.621 
19 1.1500 -0.00749 0.00749 206.579 206.723 
20 1.0000 0.00000 0.00000 234.057* 234.161* 


* Extremely close 
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Table I7. Vertical deflections, y,, of the tapered cantilever tubular beam calculated from Nastran linear 


analysis and from Shifted Displacement Transfer Functions [eq. (6b)] using the known depth factors, 
c,(=h, / 2), and Nastran linear lower surface strains, €, ; P=600 Ib at the beam tip. 


€, , in/in €,, in/in Y;5 Mi. y, ,in 
l c;, In Nastran-linear upper Nastran-linear lower Nastran-linear Shifted DTF 
surface strain surface strain vertical deflection Eq. (6b) 
0 4.0000 -0.01485 0.01485 0.000 0.000 
1 3.8500 -0.01523 0.01523 0.430 0.427 
2 3.7000 -0.01562 0.01562 1.751 1.744 
3 3.5500 -0.01603 0.01603 4.021 4.011 
4 3.4000 -0.01645 0.01645 7.307 7.296 
5 3.2500 -0.01688 0.01688 11.682 11.669 
6 3.1000 -0.01731 0.01731 17.224 17.212 
7 2.9500 -0.01775 0.01775 24.022 24.012 
8 2.8000 -0.01819 0.01819 32.174 32.166 
9 2.6500 -0.01861 0.01861 41.786 41.784 
10 2.5000 -0.01901 0.01901 52.978 52.982 
11 2.3500 -0.01937 0.01937 65.879 65.891 
12 2.2000 -0.01964 0.01964 80.632 80.656 
13 2.0500 -0.01979 0.01979 97.390 97.429 
14 1.9000 -0.01975 0.01975 116.316 116.372 
15 1.7500 -0.01940 0.01940 137.573 137.651 
16 1.6000 -0.01857 0.01857 161.313 161.415 
17 1.4500 -0.01696 0.01696 187.646 187.773 
18 1.3000 -0.01406 0.01406 216.580 216.729 
19 1.1500 -0.00899 0.00899 247.894 248 .047 
20 1.0000 0.00000 0.00000 280.868* 280.970* 
*Extremely close; beam-tip deflection = 94% of 300-in span. 
Nonlinear Analysis 


Tables [8—I14 list the strain and deflection outputs of Nastran nonlinear analysis. Nastran nonlinear 
deflection outputs have both axial and vertical displacement, {u,,y,} (fig 6b). Therefore, Nastran curved 
deflections, y, were calculated from equation (20). Note also from tables I8—I14 that the magnitudes of 


lower surface strains are slightly larger than the magnitudes of the associated upper surface strains 
because of the curved-beam effect, inducing slight axial strain components under nonlinear bending. 


In the last columns of Tables I8—114, the data of curved deflection, Y, were calculated from the 
Curved Displacement Transfer Functions [eq. (18b)] using the known depth factors, c,(= h,/2), and the 
true bending strains, (€, — €,)/2 [eq. (23)], for eliminating axial strains induced by curved-beam effect in 
nonlinear bending. 
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Table I8. Curved deflections, y,, of the tapered cantilever tubular beam calculated from Nastran 


nonlinear analysis [eq. (20)] and from Curved Displacement Transfer Functions [eq. (18b)] using the 
known depth factors, c,;(= h,/2), and the true bending strains, (€, — €,)/2 ; P=50 Ib at the beam tip. 


€, , in/in €,, in/in u,.in y,,in y, ,in y, .in 

I C; > 1 3 i> in Nesta Medios Nastran Nastran cls cneved Curved DTF 
nonlinear upper nonlinear lower : ; deflection 

surface strain surface strain x-disp. walls Eq. (20) Beeb) 
0 4.00 0.0 -0.00123 0.00123 0.000 0.000 0.000 0.000 
1 3.85 15.0 -0.00126 0.00126 0.000 0.036 0.036 0.035 
2 3.70 | 30.0 -0.00130 0.00130 0.000 0.145 0.145 0.144 
3 3.55 | 45.0 -0.00133 0.00133 0.002 0.333 0.333 0.332 
+ 3.40 | 60.0 -0.00136 0.00136 0.004 0.606 0.606 0.605 
5 3.25 | 75.0 -0.00140 0.00140 0.008 0.969 0.969 0.967 
6 3.10 | 90.0 -0.00143 0.00143 0.016 1.428 1.428 1.427 
7 2.95 | 105.0 -0.00147 0.00147 0.026 1.991 1.991 1.990 
8 2.80 |) 120.0 -0.00150 0.00151 0.041 2.666 2.666 2.666 
9 2.65 | 135.0 -0.00154 0.00154 0.062 3.462 3.463 3.462 
10 2.50 | 150.0 -0.00157 0.00157 0.091 4.388 4.389 4.390 
11 2.35 | 165.0 -0.00160 0.00160 0.129 5.455 5.457 5.458 
12 2.20 | 180.0 -0.00162 0.00162 0.179 6.675 6.679 6.680 
13 2.05 | 195.0 -0.00163 0.00163 0.243 8.059 8.064 8.068 
14 1.90 | 210.0 -0.00162 0.00163 0.324 9.621 9.629 9.634 
15 1.75 | 225.0 -0.00159 0.00160 0.427 11.374 11.386 11.393 
16 1.60 | 240.0 -0.00152 0.00153 0.555 13.330 13.347 13.356 
17 1.45 | 255.0 -0.00139 0.00140 0.712 15.497 15.521 15.532 
18 1.30 | 270.0 -0.00115 0.00116 0.902 17.875 17.909 17.921 
19 1.15 | 285.0 -0.00073 0.00074 1.124 20.445 20.491 20.504 
20 1.00 | 300.0 0.00001 0.00001 1.369 23.150 23.210* 23 .220* 


*Extremely close 
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Table 19. Curved deflections, y,, of the tapered cantilever tubular beam calculated from Nastran 


nonlinear analysis [eq. (20)] and from Curved Displacement Transfer Functions [eq. (18b)] using the 
known depth factors, c,(= h,/2), and the true bending strains, (€, — €,)/2 ; P=100 lb at the beam tip. 


é€, , in/in €, , in/in u,, in y, ,in y,,in y, ,in 
E Ci» in | x i> in Pies oe Nastran Nastran Been Sune Curved DTF 
nonlinear upper | nonlinear lower me ; deflection 
surface strain surface strain x-digp. yeaisp: Eq. (20) Eql's)) 
0 4.00 0.0 -0.00243 0.00243 0.000 0.000 0.000 0.000 
1 3.85 15.0 -0.00249 0.00249 0.000 0.070 0.070 0.070 
2 3.70 | 30.0 -0.00255 0.00255 0.002 0.287 0.287 0.285 
3 3.55 | 45.0 -0.00262 0.00262 0.006 0.658 0.658 0.656 
+ 3.40 | 60.0 -0.00268 0.00268 0.016 1.195 1.195 1.192 
5 3.25 | 75.0 -0.00275 0.00275 0.033 1.910 1.910 1.907 
6 3.10 | 90.0 -0.00281 0.00282 0.060 2.814 2.815 2.811 
7 2.95 | 105.0 | -0.00288 0.00288 0.101 3.922 3.924 3.920 
8 2.80 | 120.0 | -0.00294 0.00295 0.160 5.249 5.253 5.249 
9 2.65 | 135.0 | -0.00301 0.00301 0.242 6.812 6.818 6.815 
10 2.50 | 150.0 | -0.00306 0.00307 0.352 8.628 8.639 8.636 
11 2.35 | 165.0 | -0.00311 0.00312 0.498 10.718 10.735 10.733 
12 2.20 | 180.0 | -0.00315 0.00316 0.689 13.102 13.129 13.129 
13 2.05 | 195.0 | -0.00316 0.00317 0.934 15.803 15.844 15.847 
14 1.90 | 210.0 | -0.00314 0.00316 1.245 18.845 18.906 18.912 
15 1.75 | 225.0 | -0.00307 0.00309 1.637 22.249 22.338 22.350 
16 1.60 | 240.0 | -0.00293 0.00295 2122 26.036 26.164 26.182 
17 1.45 | 255.0 | -0.00266 0.00268 2.717 30.219 30.400 30.425 
18 1.30 | 270.0 | -0.00219 0.00222 3.431 34.793 35.043 35.077 
19 1.15 | 285.0 | -0.00139 0.00143 4.264 39.722 40.059 40.099 
20 1.00 | 300.0 0.00002 0.00002 5.184 44.896 45.337* 45 .374* 


*Extremely close 
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Table 110. Curved deflections, y,, of the tapered cantilever tubular beam calculated from Nastran 


nonlinear analysis [eq. (20)] and from Curved Displacement Transfer Functions [eq. (18b)] using the 
known depth factors, c,(= h, /2), and the true bending strains, (€, — €,)/2 ; P=200 lb at the beam tip. 


é€, , in/in €, , in/in u,, in y, ,in y, ,in j, .in 
surface strain surface strain x-disp. yedisp: Eq. (20) Bq. (heh) 
0 4.00 0.0 -0.00467 0.00467 0.000 0.000 0.000 0.000 
1 3.85 15.0 -0.00477 0.00477 0.001 0.135 0.135 0.134 
2 3.70 | 30.0 -0.00487 0.00488 0.006 0.549 0.549 0.547 
3 3.55 | 45.0 -0.00498 0.00498 0.023 1.259 1.259 1.256 
4 3.40 | 60.0 -0.00509 0.00509 0.058 2.284 2.285 2.282 
5 3.25 75.0 -0.00519 0.00520 0.120 3.644 3.647 3.644 
6 3.10 | 90.0 -0.00530 0.00531 0.219 5.362 5.369 5.366 
7 2.95 | 105.0 | -0.00540 0.00541 0.366 7.462 7.475 7AT3 
8 2.80 | 120.0 | -0.00549 0.00551 0.577 9.968 9.993 9.993 
9 2.65 | 135.0 | -0.00558 0.00560 0.868 12.909 12.953 12.955 
10 2.50 | 150.0 | -0.00565 0.00567 1.259 16.313 16.386 16.391 
11 2.35 | 165.0 | -0.00570 0.00573 1.774 20.210 20.326 20.337 
12 2.20 | 180.0 | -0.00572 0.00575 2.440 24.633 24.813 24.831 
13 2.05 | 195.0 | -0.00569 0.00574 3.290 29.611 29.883 29.910 
14 1.90 | 210.0 | -0.00561 0.00566 4.360 35.175 35.576 35.616 
15 1.75 | 225.0 | -0.00543 0.00550 5.690 41.352 41.933 41.988 
16 1.60 | 240.0 | -0.00512 0.00520 7.323 48.159 48.983 49.060 
17 1.45 | 255.0 | -0.00460 0.00469 9.298 55.600 56.748 56.853 
18 1.30 | 270.0 | -0.00375 0.00386 11.642 63.652 65.221 65.358 
19 1.15 | 285.0 | -0.00234 0.00247 14.343 72.241 74.336 74.503 
20 1.00 | 300.0 0.00008 0.00008 17.304 81.191 83.903* 84.082* 


*Extremely close 


72 


Table I11. Curved deflections, y,, of the tapered cantilever tubular beam calculated from Nastran 


nonlinear analysis [eq. (20)] and from Curved Displacement Transfer Functions [eq. (18b)] using the 
known depth factors, c,(= h, /2), and the true bending strains, (€, — €,)/2 ; P=300 Ib at the beam tip. 


é, , in/in €, , in/in u,, in y, ,in y, ,in j, ,in 
q C; > ia) i> in Nectee Se Nastran Nastran Newt cueved Curved DTF 
nonlinear upper | nonlinear lower ae : deflection 
surface strain surface strain x-disp. disp: Eq. (20) Eq (I8b) 
0 4.00 0.0 -0.00665 0.00665 0.000 0.000 0.000 0.000 
1 3.85 15.0 -0.00678 0.00678 0.001 0.193 0.193 0.191 
2 3.70 | 30.0 -0.00691 0.00691 0.013 0.781 0.781 0.778 
3 3.55 | 45.0 -0.00703 0.00704 0.047 1.789 1.790 1.785 
4 3.40 | 60.0 -0.00715 0.00716 0.117 3.241 3.244 3.238 
5 3.25 | 75.0 -0.00727 0.00728 0.241 5.164 5.172 5.166 
6 3.10 | 90.0 -0.00737 0.00740 0.437 7.585 7.604 7.596 
7 2.95 | 105.0 | -0.00747 0.00750 0.729 10.532 10.570 10.563 
8 2.80 | 120.0 | -0.00756 0.00759 1.144 14.037 14.107 14.101 
9 2.65 | 135.0 | -0.00762 0.00766 1.713 18.130 18.251 18.248 
10 2.50 | 150.0 | -0.00766 0.00771 2.472 22.843 23.043 23.044 
11 2.35 | 165.0 | -0.00766 0.00772 3.464 28.206 28.523 28.531 
12 2.20 | 180.0 | -0.00761 0.00769 4.734 34.249 34.736 34.754 
13 2.05 | 195.0 | -0.00750 0.00759 6.338 40.999 41.728 41.759 
14 1.90 | 210.0 | -0.00731 0.00742 8.333 48.476 49.539 49.591 
15 1.75 | 225.0 | -0.00699 0.00712 10.783 56.693 58.212 58.294 
16 1.60 | 240.0 | -0.00650 0.00665 13.750 65.651 67.779 67.900 
17 1.45 | 255.0 | -0.00576 0.00594 17.288 75.329 78.255 78.425 
18 1.30 | 2700 | -0.00462 0.00484 21.426 85.675 89.617 89.846 
19 1.15 | 285.0 | -0.00283 0.00309 26.133 96.589 101.780 102.065 
20 1.00 | 300.0 0.00015 0.00015 31.246 107.872 114.508* 114.820* 


*Extremely close 
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Table 112. Curved deflections, y,, of the tapered cantilever tubular beam calculated from Nastran 


nonlinear analysis [eq. (20)] and from Curved Displacement Transfer Functions [eq. (18b)] using the 
known depth factors, c,(= h, /2), and the true bending strains, (€, — €,)/2 ; P=400 lb at the beam tip. 


é€, , in/in €, , in/in u,, in y, ,in y, ,in y, in 
I C;, in Xi> in Nesta tare Nastran Nastran Neue cueve# Curved DTF 
nonlinear upper | nonlinear lower as ; deflection 
surface strain surface strain x-disp. yedisp: Eq. (20) Ea (18D) 
0 4.00 0.0 | -0.00843 0.00843 0.000 0.000 0.000 0.000 
1 3.85 15.0 | -0.00857 0.00857 0.002 0.244 0.244 0.241 
2 3.70 30.0 | -0.00870 0.00870 0.020 0.988 0.988 0.984 
3 3.55 45.0 | -0.00882 0.00883 0.074 2.260 2.262 2.255 
4 3.40 60.0 | -0.00893 0.00895 0.186 4.088 4.094 4.086 
5 3.25 75.0 | -0.00903 0.00906 0.381 6.502 6.519 6.509 
6 3.10 90.0 | -0.00911 0.00915 0.690 9.531 9.568 9.558 
7 2.95 | 105.0 | -0.00918 0.00923 1.147 13.205 13.280 13.270 
8 2.80 | 120.0 | -0.00922 0.00928 1.791 17.555 17.692 17.684 
9 2.65 | 135.0 | -0.00923 0.00930 2.668 22.610 22.845 22.841 
10 2.50 | 150.0 | -0.00919 0.00928 3.829 28.398 28.783 28.785 
11 2.35 | 165.0 | -0.00911 0.00922 5.331 34.943 35.549 35.559 
12 2.20 | 180.0 | -0.00897 0.00909 7.238 42.264 43.186 43.211 
13 2.05 | 195.0 | -0.00874 0.00888 9.619 50.376 51.740 51.785 
14 1.90 | 210.0 | -0.00841 0.00858 12.548 59.282 61.251 61.325 
15 1.75 | 225.0 | -0.00793 0.00814 16.098 68.974 71.753 71.868 
16 1.60 | 240.0 | -0.00728 0.00751 20.339 79.429 83.273 83.441 
17 1.45 | 255.0 | -0.00635 0.00663 25.326 90.600 95.812 96.048 
18 1.30 | 270.0 | -0.00502 0.00534 31.077 102.412 109.336 109.649 
19 1.15 | 285.0 | -0.00302 0.00340 37.538 114.747 123.743 124.125 
20 1.00 | 300.0 0.00022 0.00022 44.496 127.413 138.776* 139.187* 


*Extremely close 
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Table 113. Curved deflections, y,, of the tapered cantilever tubular beam calculated from Nastran 


nonlinear analysis [eq. (20)] and from Curved Displacement Transfer Functions [eq. (18b)] using the 
known depth factors, c,(= h, /2), and the true bending strains, (€, — €,)/2 ; P=500 lb at the beam tip. 


€, , in/in €, , in/in u,, in y, ,in y,,in y,.in 

E Ci» in | x i> in tere ie Nastran Nastran ne eee Curved DTF 
nonlinear upper | nonlinear lower gi : deflection 

surface strain surface strain x-disp. yodiep. Eq. (20) Eq (ts9) 
0 4.00 0.0 | -0.01005 0.01005 0.000 0.000 0.000 0.000 
1 3.85 15.0 | -0.01018 0.01019 0.003 0.291 0.291 0.288 
2 3.70 30.0 | -0.01030 0.01031 0.029 1.176 1.177 1.170 
3 3.55 45.0 | -0.01040 0.01042 0.105 2.685 2.688 2.680 
4 3.40 60.0 | -0.01049 0.01052 0.262 4.849 4.860 4.849 
5 3.25 75.0 | -0.01055 0.01059 0.534 7.698 7.726 7.714 
6 3.10 90.0 | -0.01059 0.01065 0.963 11.261 11.323 11.310 
7 2.95 | 105.0 | -0.01060 0.01067 1.594 15.568 15.690 15.678 
8 2.80 | 120.0 | -0.01058 0.01066 2.479 20.644 20.867 20.856 
9 2.65 | 135.0 | -0.01051 0.01061 3.674 26.514 26.894 26.887 
10 2.50 | 150.0 | -0.01038 0.01051 5.245 33.198 33.816 33.815 
11 2.35 | 165.0 | -0.01019 0.01034 7.259 40.708 41.672 41.682 
12 2.20 | 180.0 | -0.00993 0.01010 9.792 49.052 50.505 50.532 
13 2.05 | 195.0 | -0.00957 0.00977 12.922 58.225 60.353 60.404 
14 1.90 | 210.0 | -0.00909 0.00933 16.729 68 .213 71.251 71.336 
15 1.75 | 225.0 | -0.00847 0.00875 21.288 78.986 83.225 83.356 
16 1.60 | 240.0 | -0.00765 0.00798 26.667 90.496 96.286 96.480 
17 1.45 | 255.0 | -0.00658 0.00696 32.909 102.677 110.428 110.695 
18 1.30 | 270.0 | -0.00512 0.00556 40.018 115.437 125.603 125.950 
19 1.15 | 285.0 | -0.00302 0.00353 47.913 128.652 141.694 142.107 
20 1.00 | 300.0 0.00030 0.00030 56.353 142.146 158.439* 158.866* 


*Extremely close 
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Table 114. Curved deflections, y,, of the tapered cantilever tubular beam calculated from Nastran 


nonlinear analysis [eq. (20)] and from Curved Displacement Transfer Functions [eq. (18b)] using the 
known depth factors, c,(= h,/2), and the true bending strains, (€, — €,)/2 ; P=600 lb at beam tip. 


€, , in/in €,, in/in u,, in y, ,in y, ,in 3, ,in 

I Ci» in | x i> in soins haaees Nastran Nastran es paves Curved DTF 
nonlinear upper | nonlinear lower : : deflection 

surface strain surface strain aulsp yep. Eq. (20) Bae) 
0 4.00 0.0 | -0.01154 0.01155 0.000 0.000 0.000 0.000 
1 3.85 15.0 | -0.01166 0.01167 0.004 0.334 0.334 0.330 
2 3.70 30.0 | -0.01175 0.01177 0.038 1.348 1.349 1.342 
3 3.55 45.0 | -0.01182 0.01185 0.137 3.074 3.079 3.069 
4 3.40 60.0 | -0.01187 0.01191 0.342 5.543 5.559 5.545 
5 3.25 75.0 | -0.01189 0.01195 0.695 8.783 8.824 8.809 
6 3.10 90.0 | -0.01187 0.01194 1.249 12.824 12.915 12.898 
7 2.95 | 105.0 | -0.01181 0.01190 2.059 17.689 17.868 17.851 
8 2.80 | 120.0 | -0.01170 0.01182 3.188 23.399 23.724 23.708 
9 2.65 | 135.0) -0.01154 0.01168 4.703 29.971 30.522 30.509 
10 2.50 | 150.0 | -0.01131 0.01148 6.678 37.412 38.301 38.296 
11 2.35 | 165.0 | -0.01100 0.01120 9.191 45.725 47.100 47.107 
12 2.20 | 180.00 | -0.01061 0.01084 12.322 54.901 56.955 56.980 
13 2.05 | 195.0 | -0.01011 0.01038 16.154 64.917 67.897 67.949 
14 1.90 | 210.0 | -0.00949 0.00980 20.765 75.740 79.950 80.040 
15 1.75 | 225.0 | -0.00873 0.00909 26.227 87.320 93.131 93.271 
16 1.60 | 240.0 | -0.00778 0.00819 32.597 99.593 107.442 107.643 
17 1.45 | 255.0 | -0.00658 0.00706 39.905 112.474 122.860 123.132 
18 1.30 | 270.0 | -0.00504 0.00559 48.133 125.863 139.326 139.668 
19 1.15 | 285.0 | -0.00291 0.00355 57.183 139.637 156.715 157.104 
20 1.00 | 300.0 0.00037 0.00037 66.792 153.640 174.763* 175.138* 


* 


Extremely close; beam-tip deflection = 58% of 300-in span. 
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